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Chapter  1.  Introduction 


This  report  is  concerned  with  propagation,  between  the  Earth's 
surface  and  the  ionosphere,  of  electromagnetic  waves  having  frequencies 
much  lower  than  those  customarily  used  for  radio  communications.  The 
calculations  refer  to  the  frequency  range  7.5-300  Hz.  At  these  fre¬ 
quencies  the  height  of  the  ionosphere  is  much  less  than  the  wavelength 
in  free  space.  Consequently  the  Earth  and  the  ionosphere  constitute 
the  two  conductors  of  a  transmission  line.  The  problem  is  to  calculate 
the  effect  of  ionospheric  phenomena  on  the  propagation  properties  of 
such  a  transmission  line,  and  to  understand  the  physical  reasons  for  the 
results  obtained. 

Use  of  extremely  low  frequencies  for  practical  radio  communica¬ 
tions  is  discussed  in  a  special  issue  of  the  IEEE  Transactions  on 
Communications  edited  by  Wait  (1974).  Such  a  communications  system 
would  be  able  to  penetrate  considerably  further  into  the  ocean  than  is 
possible  at  conventional  radio  frequencies. 

Another  reason  for  interest  in  extremely  low  radio  frequencies 
is  that  the  space  between  the  conducting  Earth  and  the  conducting 
ionosphere  is  a  cavity  resonator  of  enormous  dimensions.  The  funda¬ 
mental  wavelength  is  controlled  by  the  perimeter  of  the  Earth,  and 
corresponds  to  a  frequency. of  about  8  Hz.  The  existence  of  this 
resonance  was  first  pointed  out  by  Schumann  (1952a, b,c).  It  is 
excited  by  electromagnetic  radiation  from  lightning  flashes,  dis¬ 
tributed  over  the  world.  Such  radiation  constitutes  atmospheric 
noise,  and  there  are  peaks  of  noise  at  the  Schumann  resonance 
frequency  and  at  its  low  harmonics  (Schumann  and  Konig  1954, 
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Schumann  1957,  Balser  and  Wagner  1962a, b,  Chapman  and  Jones  1964, 

Madden  and  Thompson  1965,  Jones  1967,  Larsen  and  Egeland  1968).  The 
strength  of  these  noise  peaks,  the  precise  frequencies  at  which  they 
occur,  and  their  bandwidths  are  related  to  the  properties  of  the 
ionosphere.  It  is  desirable  to  understand  these  features,  and  it  may 
even  be  possible  to  use  the  structure  of  the  Schumann  resonances  to 
derive  some  information  about  the  physics  of  the  ionosphere. 

For  many  years  the  theory  of  ELF  propagation  in  the  Earth- 
ionosphere  transmission  line  was  based  on  the  assumption  that  the 
ionosphere  may  be  regarded  as  one  or  more  homogeneous  conducting  slabs 
with  discontinuous  boundaries  at  heights  chosen  in  a  manner  that  was 
somewhat  arbitrary  but  nevertheless  reasonable  (Watson  1919,  Budden 
1961a ,b.  Wait  1962,  Galejs  1972).  The  scale-height  of  the  ionosphere  is 
small  compared  with  the  free-space  wavelength,  and  in  this  sense  the 
under  boundary  of  the  ionosphere  is  sharp.  Nevertheless,  a  need 
existed  to  employ  a  continuous  vertical  profile  of  conductivity  that 
had  some  relation  to  the  known  facts  of  ionospheric  physics.  A  need 
also  existed  to  take  full  account  of  the  effect  of  the  Earth's 
magnetic  field. 

Because  of  the  long  free-space  wavelength,  the  phase- integral 
treatment  (WKB  approximation)  cannot  be  used  near  the  bottom  of  the 
ionosphere.  On  the  other  hand,  full-wave  treatments  (Pappert  and 
Moler  1974)  lead  to  elaborate  computer  programs.  These  give  numerical 
values  for  the  propagation  parameters  in  the  Earth-ionosphere  trans¬ 
mission  line,  but  only  limited  insight  into  the  physics  of  the 
propagation  phenomena  in  the  ionosphere  (see  Pappert  and  Moler  1978) . 
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»A.  » 


A  need  exists  for  an  approximate  treatment  of  ELF  propagation 
that  reveals  the  physics  of  what  is  happening  in  the  ionosphere. 

Such  a  theory  needs  to  take  into  account  the  ionospheric  profiles  of 
the  electron  density,  of  the  ion  densities,  and  of  the  collisional 
frequencies,  as  well  as  the  effect  of  the  Earth's  magnetic  field.  It 
needs  to  be  such  that  a  particular  propagation  phenomenon  can  be 
traced  to  the  feature  of  the  ionospheric  profiles  that  is  principally 
responsible.  Such  a  theory  was  attempted  several  years  ago  by 
Booker  and  Lefeuvre  (1977) ,  and  their  method  is  described  in  Chapter  2 
of  this  report.  The  theory  shows  that  ELF  waves  are  reflected  in  a 
stratum  where  1/(4tt)  times  the  local  wavelength  is  of  the  same  order  of 
magnitude  as  the  local  scale-height  of  the  ionospheric  profile. 

However,  the  theory  of  Booker  and  Lefeuvre  had  two  drawbacks. 

Firstly,  it  did  not  take  satisfactory  account  of  the  effect  of 
ionization  below  the  reflecting  stratum.  A  method  for  doing  this  has 
since  been  developed  by  Greifinger  and  Greifinger  (1978,  1979). 
Secondly,  the  theory  of  Booker  and  Lefeuvre  assumed  that  only  one 
reflecting  stratum  exists  in  the  ionosphere  whereas  several  are 
frequently  encountered.  By  day,  reflection  occurs  both  from  the 
gradient  on  the  under  side  of  the  D  region  and  from  the  gradient  on 
the  under  side  of  the  E  region.  By  night,  reflection  occurs  from 
the  gradient  on  the  top  side  of  the  E  region  as  well  as  from  the 
gradient  on  the  under  side  of  the  E  region.  Below  the  ionic  gyro- 
frequency  by  night,  there  is  even  reflection  from  the  F  region. 

The  theory  of  Booker  and  Lefeuvre  (1977)  has  been  extended 
in  order  (i)  to  cover  the  occurence  of  several  reflecting  strata, 
and  (ii)  to  incorporate  the  Greifinger  and  Greifinger  (1978,  1979) 
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treatment  of  ionization  near  the  bottom  of  the  ionosphere.  The 
extended  theory  is  described  in  Chapter  3  of  this  report. 


On  the  basis  of  this  theory,  calculations  have  been  carried  out 
by  Amir  B.  Behroozi-Toosi,  a  graduate  student  studying  for  a  Master's 
degree.  A  simplified  worldwide  model  of  the  ionosphere  is  employed 
to  describe  approximately  the  behavior  of  the  D  and  E  regions  as 
functions  of  latitude,  time  of  day,  season  of  the  year,  and  year  in  the 
sunspot  cycle.  A  central  dipole  model  is  used  for  the  Earth's  magnetic 
field.  The  F  region  was  omitted  on  fhe  assumption  that  it  is  not 
important  for  ELF  propagation  in  the  Earth- ionosphere  transmission 
line — an  assumption  that  proved  to  be  less  true  than  anticipated.  On 
the  basis  of  the  theory  presented  in  Chapter  3,  calculations  were 
made  for  the  frequency  range  7.5-300  Hz,  and  results  are  presented  in 
the  diagrams  of  Chapter  4.  These  are  designed  so  as  not  only  to 
describe  the  propagation  characteristics  of  the  Earth- ionosphere 
transmission  line  but  also  to  bring  out  the  physical  reasons  for  the 
results  obtained. 

A  comparison  is  made  between  the  results  obtained  by  the  method 
of  Chapter  3  and  those  obtained  by  the  method  of  Pappert  and  Moler 
(1974) .  The  comparison  is  presented  in  Tables  4  and  5  on  page  153, 
the  latter  table  being  provided  through  the  courtesy  of  the  Naval 
Ocean  Systems  Center.  At  present,  and  probably  for  a  long  time  to 
come,  overall  error  in  such  calculations  will  be  controlled  by 
ignorance  concerning  worldwide  behavior  of  the  ionosphere.  For  this 
reason,  little  is  gained  by  going  from  moderate  accuracy  to  high 
accuracy  in  the  propagation  calculations.  On  the  other  hand,  much  is 
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gained  by  understanding  what  features  of  the  ionospheric  model  control 
what  propagation  phenomena.  It  is  believed  that  the  method  described 
in  Chapter  3  and  applied  in  Chapter  4  provides  just  what  is  required 
It  gives  more  physical  insight  than  any  other  method  currently  avail¬ 
able.  It  uses  less  effort  than  any  other  profile-based  method  cur¬ 
rently  available.  In  spite  of  approximations,  it  provides  reasonable 
accuracy . 
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Chapter  2.  The  reflection  theory  of 
Booker  and  Lefeuvre  (1977) 


1.  Introduction 

Propagation  of  the  transmission-line  mode  in  the  wave-guide 
between  the  conducting  Earth  and  the  conducting  ionosphere  was  first 
studied  by  Watson  (1919),  and  the  theory  has  since  been  developed,  par¬ 
ticularly  for  application  in  the  ELF  band.  It  is  summarized  in  books 
by  Budden  (1961b),  Wait  (1962)  and  Galejs  (1972),  and  also  in  a  special 
issue  of  the  IEEE  Transactions  on  Communications  edited  by  Wait  (1974) . 
The  earth,  the  atmosphere  and  the  ionosphere  are  often  regarded  as  three 
media  that  are  homogeneous  and  isotropic.  More  elaborate  versions  of 
the  theory  take  the  ionosphere  to  involve  two  or  more  homogeneous  iso¬ 
tropic  layers  (Wait,  1958;  Chapman  and  Jones,  1964a;  Jones,  1967). 
Theories  involving  a  more-or-less  continuously  stratified  ionosphere, 
isotropic  or  non-isotropic,  have  been  developed  mainly  for  VLF  and 
higher  frequencies,  but  they  have  also  been  used  at  ELF  (Cole,  1965; 
Madden  and  Thompson,  1965;  Hughes  and  Theisen,  1970;  Barr,  1972; 

Pappert  and  Moler,  1974). 

Even  if  it  is  agreed  that  a  homogeneous  ionosphere  with  a  dis¬ 
continuous  under  boundary  is  an  adequate  model  at  ELF,  we  need  to  know 
at  what  height  to  place  the  discontinuity  and  what  parameters  to  assign 
to  the  ionosphere  above  this  height.  Furthermore,  from  experience  at 
higher  frequencies,  one  would  expect  the  height  of  reflection  and  the 
other  relevant  parameters  to  depend  on  angle  of  incidence  and  frequency. 
Konrad  and  Poverlein  (1958)  and  Chapman  and  Jones  (1964b)  have  con¬ 
jectured  that  the  level  of  reflection  should  be  higher  at  ELF  than  at 
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VLF.  Nevertheless,  it  is  frequently  assumed  that  the  same  ipnospheric 
discontinuity  is  satisfactory  for  all  angles  of  incidence  and  all  fre¬ 
quencies  in  the  ELF  band.  If  so,  then  it  should  be  possible  to  give 
theoretical  reasons  for  this  result.  At  the  Schumann  resonant  fre¬ 
quencies,  this  matter  has  been  considered  by  Large  and  Wait  (1968). 

It  is  the  objective  of  this  chapter  to  start  with  a  simple  but 
reasonably  realistic  profile  of  ionization  density,  and  to  derive  a 
method,  suitable  for  use  at  ELF,  for  ascertaining  at  what  level  the 
ionosphere  may  best  be  cut  off  at  the  bottom  and  be  replaced  by  a  dis¬ 
continuity.  We  shall  also  ascertain  the  dependence  of  this  discontin¬ 
uity-level  on  angle  of  incidence  and  frequency,  and  go  on  to  assess  the 
rate  of  attenuation  in  the  Earth-ionosphere  transmission  line  due  to 
leakage  of  energy  into  the  ionosphere  above  the  level  of  reflexion. 

We  shall  do  this  taking  into  account  the  effect  of  the  Earth's 
magnetic  field,  but  neglecting  coupling  between  the  ordinary  and  extra¬ 
ordinary  waves  (the  0  and  X  waves).  It  will  be  necessary  to  solve  the 
mode  transmission-problem  in  the  Earth-ionosphere  transmission  line 
taking  the  ionospheric  conductor  as  a  non-isotropic  medium.  It  should 
be  noted  that,  in  such  circumstances,  a  mode  involves,  not  two  crossing 
plane  waves,  but  four.  There  are  a  pair  of  upgoing  waves  and  a  pair  of 
downgoing  waves,  the  members  of  each  pair  having  the  elliptic  polariza¬ 
tion  corresponding  to  the  0  and  X  waves  in  the  ionosphere. 
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>.  The  level  of  reflection  in  the  ionosphere  at  ELF 

We  shall  base  our  discussion  of  the  level  of  reflection  in  the 
ionosphere  at  ELF  on  an  ionospheric  model  for  which  the  index  of  re¬ 
fraction  n  varies  with  height  2  according  to  the  formula 

n^  »  1  -  (bexpjij;)  exp  {(z-a)/h)  (1) 

where  b  and  ii  are  real  constants,  h  is  a  scale-height,  and  a  is  a 
reference  height.  Such  a  profile  of  refractive  index  is  not  ideal,  but 
gives  a  reasonable  representation  of  the  underside  of  the  ionosphere  in 
circumstances  when  the  increase  of  ionization  density  with  height  is 
roughly  exponential.  The  profile  does  not  represent  properly  the  vari¬ 
ation  of  collisional  frequency  with  height,  as  a  result  of  which  there 
is  a  tendency  to  underestimate  absorption  on  one  side  of  the  level  of 
reflection  and  to  overestimate  it  on  the  other  side. 

For  a  plane  wave  incident  upon  a  plane  stratified  ionosphere 
whose  profile  is  described  by  Equation  (1)  there  is  an  exact  wave 
solution  for  horizontal  polarization.  We  shall  assume  that  this  solu¬ 
tion  gives  a  reasonable  clue  to  determination  of  the  level  of  reflection 
in  a  practical  ionosphere  for  the  0  and  X  waves.  The  solution  is  given 
by  Budden  (1961a,  pp.  354-357)  and  by  Wait  (1962,  pp.  72-74).  It 
involves  a  Bessel  function  of  an  argument  that  is  small  at  the  bottom 
of  the  ionosphere  and  large  well  inside  the  ionosphere,  with  reflection 
occurring  in  the  transition  region  between  the  two.  For  an  angle  of 
incidence  9  (S*  sin  9,  O  cos  8),  possibly  complex,  the  change  in  com¬ 
plex  phase  between  the  incident  wave  (measured  at  the  reference  height 
a  and  the  reflected  wave  (also  measured  at  z  =  a)  is 
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where  A^  is  the  wavelength  in  free  space.  When  2Trh/\o  is  large. 

Equation  (2)  corresponds  to  the  phase-integral  (or  WKB)  approximation 

(Budden,  1961a).  On  the  other  hand,  at  ELF  2rrh/A  is  small,  and  we 
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Then  Equation  (2)  may  be  written  at  ELF 
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or 
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For  the  ionospheric  profile  in  Equation  (1)  the  phase-integral 
approximation  always  applies  at  sufficiently  large  heights.  A  suffic¬ 
iently  large  height  is  one  for  which  the  argument  of  the  Bessel  wave- 
function  is  appreciably  greater  than  unity.  Any  height  above  the  zone 
of  transition  from  small  argument  to  large  argument  satisfies  this 
condition,  and  this  zone  is  the  reflecting  stratum.  Consequently  we 
may  simply  say  that  the  phase-integral  approximation  applies  above  the 
reflecting  stratum. 

Let  us  compare  the  result  in  Equation  (6)  with  that  for  an 
ionosphere  in  which  (i)  the  same  refractive  index  profile  is  maintained 
above  a  certain  level  z^,  (ii)  the  phase- integral  approximation  is 
applied  for  z  >  z^  and  (iii)  the  refractive  index  is  taken  as  unity  for 
z  <  z^.  Let  us  see  whether  it  is  possible  to  adjust  the  level 
z ^  of  the  resulting  discontinuity  in  refractive  index  so  that,  well 
below  the  level  z^,  the  phase-change  1  between  an  incident  wave 
(measured  at  the  reference  height  a)  and  the  reflected  wave  (measured 


also  at  2  =  a)  is  nearly  the  same  as  the  right  hand  side  of 
Equation  (6) . 

Just  above  the  level  z^  the  square  of  the  refractive  index  is 


n 


2 

1 


1 


(bexpjtjj)  exp  {(z^  -a)/h] 
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and  if  we  assume  that  this  level  is  sufficiently  far  into  the  iono¬ 
sphere  to  make  |n|»  1,  we  may  write  approximately 


n^  =  -(bexpjv)  exp  {(z^  -a)/h} 


(8) 


For  an  angle  of  incidence  whose  cosine  in  C,  the  Fresnel  reflection  co¬ 
efficient  at  the  discontinuity  at  level  z^  is 
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Assuming  that  |n^]>>l»C,  this  becomes 
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which  is  equivalent  to  a  complex  phase-change  on  reflection  of 
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and  does  not  differ  much  from  it.  This  is  the  phase-change  referred  to 
as  reference  height.  If  a  is  used  as  the  reference  height,  the  com¬ 
plex  phase-change  between  the  incident  and  reflected  waves  is 
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To  put  this  result  in  approximately  the  same  form  as  Equation  (6)  , 
we  note  that  Equation  (8)  may  be  rewritten 
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so  that  Equation  (12)  becomes 
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or,  again  using  Equation  (8), 


$  '  -  TT  + 


4tt  hC 

X 


In 


2C  /  -J  4  \ 

— ,  (cosrJS1nT) 


(14) 


We  are  now  in  a  position  to  deduce  from  Equations  (6)  and  (14) 
whether  it  is  possible  to  choose  z^  (i.e.  n^)  so  as  to  make  the  complex 
phase-change  i’  for  the  discontinuity  profile  approximately  the  same  as 
that  for  the  original  profile.  We  see  that  it  is  not  possible  to  do  so 
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precisely  because  of  the  terms  following  the  logarithms  in  Equations 

(6)  and  (14).  On  the  other  hand,  we  also  see  that  these  terms  are 

small  compared  with  the  logarithmic  terms  when  A  / (2Th)  and  In,  j  are 

o  1 

large,  that  is,  at  ELF.  To  make  the  dominant  logarithmic  terms  equal 
we  require  that 


A 


This  means  that  the  reflecting  stratum  in  the  ELF  band  occurs 
close  to  where  the  local  reduced  wavelength  A  /(2T|n.j)  is  equal  to  the 
scale  height  h  of  the  profile.  Appreciably  below  the  level  correspond¬ 
ing  to  Equation  (15)  the  local  wavelength  is  large  compared  with  the 
structure-scale  of  the  ionization  prof ile,  which  is  therefore  largely 
ignored  by  the  waves.  Appreciably  above  the  level  corresponding  to 
Equation  (15)  the  local  wavelength  is  small  compared  with  the  structure- 
scale  of  the  ionization  profile  and  the  phase-integral  approximation  is 
applicable.  It  is  in  the  neighborhood  of  the  level  given  by  Equation 
(15)  that  the  bulk  of  the  reflection  takes  place. 

Let  us  now  examine  the  small  terms  following  the  logarithms  in 
Equations  (6)  and  (14) .  While  it  is  not  in  general  possible  to  adjust 
the  level  z^  so  as  to  make  these  terms  exactly  equal,  it  is  possible  to 
adjust  the  level  in  such  as  way  as  to  minimize  the  difference  between 
P  and  .  From  Equations  (6)  and  (14)  we  deduce  that 
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Figure  1.  Illustrating  (i)  the  value  of  y  required  to  minimize  error 
in  the  complex  phase-change  on  reflection,  and  (ii)  the 
value  of  S  to  be  used  in  expression  (19)  for  estimating 
the  resulting  error. 
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For  a  value  of  y  differing  somewhat  from  unity,  Equation  (18)  deter¬ 
mines  a  height  a  little  different  from  that  given  by  Equation  (15).  We 
select  y  so  as  to  minimize  the  expression  for  j  ■*>  —  <t> * )  *"  given  in 
Equation  (17).  This  value  of  y  turns  out  to  be  independent  of  fre¬ 
quency,  of  scale-height  and  of  angle  incidence,  but  to  be  dependent  on 
the  phase-angle  p  of  the  refractive  index  profile  in  Equation  (1). 

The  value  of  y  varies  between  the  extreme  values  2.00  and  2.65  as  shown 
in  Figure  1,  and  even  less  in  practice.  With  the  value  of  y  chosen  in 
this  way,  the  error  in  |$  -  •*> ’  J  has  magnitude 


18 


(19) 
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where  6  is  as  shown  in  Figure  1.  This  expression  has  Co  be  compared 
with  Che  logarithmic  term  in  Equations  (6)  or  (14) ,  and  therefore  with 
the  term 


(z1  -  a)  (20) 

'  o 

in  Equation  (12).  The  condition  that  4>  may  be  replaced  by  ’  in  a 

mode  calculation  involving  ionospheric  reflexion  at  a  height  z  above 

K 

the  surface  of  the  earth  is  therefore 


4T:hC 

A 

o 


S  « 


4ttC 

A 


o 


z 


R 


or 

3  «  zR/h  (21) 

Since  z_/h  is  of  the  order  of  20  or  more,  and  the  maximum  value  of  g 
K 

in  Figure  1  is  1.77,  it  is  clear  that  the  condition  (21)  is  reasonably 
well  satisfied. 

This  implies  that,  to  a  first  approximation,  one  may  abolish 
the  ionosphere  below  the  level  where 
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2trh 


(22) 


and  apply  the  phase-integral  approximation  above  this  level.  The 
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quantity  y  is  given  by  Figure  1  and  is  independent  of  the  angle  of 
incidence.  It  follows  that  the  level  at  which  a  discontinuous  iono¬ 


spheric  lower-boundary  should  be  placed  for  ELF  studies  is  in  fact  in¬ 
dependent  of  the  angle  of  incidence  (including  complex  angles  of  in¬ 
cidence).  On  the  other  hand,  the  discontinuity-level  does  depend  on 
frequency:  not  only  does  Xq  appear  explicitly  in  Equation  (22),  but 
the  refractive  index  n  is  also  frequency  dependent.  Hence  it  is  not 
advisable  to  use  the  same  discontinuity-level  throughout  the  ELF  band, 
although  at  each  frequency  it  is  satisfactory  to  use  the  same  discon¬ 
tinuity-level  for  all  angles  of  incidence. 

While  the  above  deductions  apply  strictly  to  the  refractive 
index  profile  in  Equation  (1),  it  is  not  difficult  to  devise  a  way  to 
apply  the  results  to  more  practical  ionospheric  profiles.  For  such  a 
profile  a  scale-height  can  be  defined  at  each  height  for  each  fre¬ 
quency  by  means  of  the  equation 


h 


1 


|-1 

1) 


(23) 


and  it  is  then  simple  to  ascertain  the  height  where  this  is  equal  to  the 
local  reduced  wavelength  J  27rn  j .  This  height  is  roughly  the  ELF 

height  of  reflection  in  accordance  with  Equation  (15).  At  this  height 

one  can  ascertain  the  value  of  the  complex  refractive  index  n^  for  the 
given  profile  and  hence  the  quantity 

P  =  tt  +  2  arg  n^  (24) 

In  accordance  with  Equation  (8) ,  this  is  approximately  the  value  of  ’p 
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to  be  used  in  Figure  1  in  order  to  deduce  a  value  for  y.  We  can  now 
ascertain  the  level  where  the  scale-height  given  by  Equation  (23)  is 
equal  to  the  corrected  local  reduced  wavelength  Aq/ j  2tt  y  n|.  In  ac¬ 
cordance  with  Equation  (22)  this  height  is  a  reasonable  approximation 
to  the  ELF  height  of  reflection  for  the  selected  frequency  and  for  any 
angle  of  incidence.  Substitution  from  Equation  (23)  into  Equation  (22) 
shows  that  the  ELF  height  of  reflection,  defined  in  this  way,  is  given 
approximately  by 


I  dn^  _  2tt  y 

3  dz  X 

n  o 


(25) 


In  the  presence  of  the  Earth's  magnetic  field,  the  above  tech¬ 
nique  may  be  applied  to  the  refractive  index  profiles  for  0  and  X 
waves,  provided  that  one  neglects  coupling  between  these"  waves.  The 
refractive  indices  at  the  levels  given  by  Equation  (25)  always  turn  out 
to  be  large  compared  with  unity.  Hence  the  ELF  waves  transmitted  into 
the  ionosphere  may  be  assumed  to  travel  almost  vertically  upwards  re¬ 
gardless  of  the  direction  of  the  incident  and  reflected  waves  below 
the  ionosphere.  This  greatly  simplifies  the  handling  of  the  Earth's 
magnetic  field.  The  refractive  indices  to  be  used  in  Equation  (25) 
for  the  0  and  X  waves  are  simply  those  for  vertical  propagation  in  the 
ionosphere. 


21 


3.  The  mode  theory  at  ELF 

We  derive  the  mode  theory  in  the  Earth- ionosphere  wave-guide 
at  ELF  on  the  assumption  that  curvature  of  the  Earth  may  be  neglected, 
that  the  height  of  reflection  is  given  by  Equation  (22),  that  the  phase- 
integral  approximation  may  be  applied  above  this  level,  and  that  the 
ionosphere  may  be  abolished  below  this  level.  Reflection  at  the  result¬ 
ing  discontinuity  is  then  determined  by  the  refractive  index  immedi¬ 
ately  above  the  discontinuity. 

We  derive  separate  heights  of  reflection  zQ  and  z for  the  0  and 
X  waves  from  Equation  (22)  using  the  refractive  indices  for  the  0  and  X 
waves  for  vertical  propagation.  Just  above  the  associated  discontinu¬ 
ities,  let  nQ  and  nx  be  the  refractive  indices  for  the  0  and  X  waves, 
so  that  the  impedances  looking  vertically  upwards  at  the  location  of  the 
discontinuities  are  C  /nn  and  t  / nv,  where  C  is  the  impedance  of  free 
space  (377  ohms).  If  we  take  axes  as  shown  in  Figure  2,  with  z  ver¬ 
tically  upwards  and  the  yz  plane  as  the  vertical  plane  of  propagation 

of  the  mode,  then  we  have,  in  terms  of  the  components  of  the  complex 

0  0  —X  X 

electromagnetic  fields  12  ,  H  for  the  0  wave  and  £  ,  H*  for  the  X  wave 
at  their  respective  levels  of  reflection. 


(26) 


(27) 
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Also,  If  Qq  and  Qx  are  the  polarization  ratios  for  the  0  and  X  waves 
at  their  respective  levels  of  reflection,  we  have 


(28) 


(29) 


A  mode  in  the  Earth-ionosphere  wave-guide  consists  of  a  pair  of 
crossing  plane  waves  at  an  angle  0  to  the  vertical  (S  =  sin  0, 

C  =  cos  6)  having  an  elliptic  polarization  whose  horizontal  projection 
corresponds  to  Equation  (28)  combined  with  a  pair  of  crossing  plane 
waves  at  the  same  angle  9  having  an  elliptic  polarization  whose  hori¬ 
zontal  projection  corresponds  to  Equation  (29).  From  the  vertical 
standing-wave  pattern  formed  by  these  waves  we  may  write  down,  by 
standard  matrix  methods  for  one-dimensional  transmission  systems  (see, 
for  example,  Booker,  1948),  the  complex  electromagnetic  field  (^,  )  at 

the  ground  in  terms  of  the  electromagnetic  fields  at  the  levels  of 
reflection.  We  obtain: 
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,  2tt  z  C  Zttz  C 

rT  1  O  ,  O  ttO 
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To  satisfy  the  mode  condition  these  values  of  E  ,  E  ,  H  and  H  must 

x  y  x  y 

obey  the  boundary  conditions  at  the  surface  of  the  earth. 

By  employing  the  ratios  (26),  (27),  (28)  and  (29)  one  can  express 
all  the  electromagnetic  field  components  on  the  right  hand  sides  of 
Equations  (30),  (31),  (32)  and  (33)  in  terms  of  and  H^,  which  are 
the  components  of  the  magnetic  field  for  the  0  and  X  waves  at  their 
levels  of  reflection,  taken  in  the  direction  transverse  to  that  of 
propagation  in  the  Earth-ionosphere  wave-guide.  We  obtain: 
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and 
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are  sufficiently  small  that  their  sines  may  be  replaced  by  the  angles 
and  their  cosines  by  unity.  Equations  ( 34 )  -  (37)  may  then  be  written 
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(40) 
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In  these  equations  the  quantities 
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(42) 


are  all  small  in  practice.  Care  is  therefore  necessary  in  making 
further  approximations  in  Equations  (38)  and  (40).  However,  Equations 
(39)  and  (41)  can  safely  be  simplified  to 
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These  equations  relate  the  two  Cartesian  components  (H^,  H  )  of  the 

horizontal  magnetic  field  of  the  wave  at  the  Earth's  surface  to  the 
0  X 

quantities  and  that  measure  the  strengths  of  the  0  and  X  waves 
at  their  respective  levels  of  reflection. 

We  now  have  to  satisfy  the  boundary  conditions  at  the  Earth's 
surface.  There  is,  of  course,  some  leakage  of  energy  into  the  Earth 
from  the  Earth- ionosphere  transmission  line.  We  know  how  to  calculate 
the  corresponding  contribution  to  attenuation  in  the  transmission  line 
(Wait  1962)  and  we  know  that  it  is  additive  to  the  ionospheric  con¬ 
tribution.  Since  we  are  at  present  interested  in  calculating  the 
ionospheric  contribution,  we  can  take  the  Earth  as  a  perfect  conductor 
and  write  the  boundary  condition  at  the  Earth's  surface  as 


E  =  E  =  0,  z  =  0 
x  y 


It  follows  from  Equations  (38)  and  (40)  that 


—  +  j 

no 


QH°+  (  -=•  +  j 


-X  \q  hx  =  0 


0  X 

Elimination  of  the  ratio  between  this  pair  of  equations  gives 


the  equation  for  the  value  of  C  appropriate  to  the  transmission-line 
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The  complex  propagation  constant  along  the  Earth-ionosphere  trans¬ 
mission  line  is  then 


where 


2rr  S 
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S  =  (1  -  C 


2.1/2 


(50) 


(51) 


and  C  is  given  by  Equation  (49 )  .  The  real  part  of  S  is  the  ratio  of 
velocity  of  light  to  the  velocity  of  phase  propagation  in  the  Earth- 
ionosphere  transmission  line,  while  the  imaginary  part  of  expression 
(50)  with  the  sign  reversed  is  the  rate  of  attenuation  in  nepers  per 
unit  distance. 

If  one  is  interested  in  the  relative  degree  of  excitation  of  the 
0  and  X  waves,  this  is  given  by  Equation  (47)  in  the  form 
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If  one  is  interested  in  the  slight  degree  of  horizontal  elliptic 
polarization  existing  at  ground  level,  one  divides  Equation  (44)  by 
Equation  (43)  and  uses  Equation  (52)  to  obtain 
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4.  Useful  analytical  approximations 


While  for  numerical  purposes  there  is  no  problem  about  making 
mode  calculations  using  Equations  (49)  -  (53),  it  should  be  noted  that, 
for  analytical  purposes,  some  further  approximations  can  usefully  be 
made.  Although  all  the  quantities  (42)  are  small,  they  are  not  of 
the  same  order  of  smallness.  In  practice  we  have 
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but  we  also  have 
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In  consequence  the  brackets  in  Equation  (40)  simplify,  whereas  those 
in  Equation  (38)  do  not,  and  a  corresponding  statement  applies  for 
Equations  (46)  and  (47).  Hence  Equation  (49)  simplifies  to 
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Moreover  Equation  (51)  may  usually  be  written  approximately  as 
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with  the  result  that  the  value  of  S  to  be  used  in  the  propagation 
constant  (50)  is 
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Likewise  Equation  (52)  for  the  relative  excitation  of  the  0  and  X 
waves  becomes 
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while  Equation  (53)  for  the  horizontal  elliptic  polarization  at  ground 
level  becomes 
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The  slight  ellipticity  in  the  horizontal  field  at  ground  level  depends 
primarily  on  the  small  difference  in  the  heights  of  reflection  for 
the  0  and  X  waves. 


5.  Cold  magneto-plasma  relations 


To  apply  the  methods  developed  above,  it  is  necessary  to  use 
expressions  for  the  complex  refractive  indices  and  polarization 
ratios  of  the  0  and  X  waves.  At  ELF,  it  is  desirable  to  take  into 
account  ions,  positive  and  negative,  atomic  and  molecular.  Five 
species  of  charge  carriers  are  therefore  involved:  electrons  (e) , 
positive  atomic  ions  (a  +) ,  negative  atomic  ions  (a  -) ,  positive  molec¬ 
ular  ions  (m  +) ,  and  negative  molecular  ions  (m  -) .  For  each  of  these 
species  we  can  define  angular  magnetic  gyro  frequencies  -!^e»  +‘i^a+» 
-ui.  ,  +lJ^jra+  and  noCice  that  the  definitions  of  u^e>  coMa_  and 

^Mm-  are  suc*1  that  they  are  numerically  negative.  For  each  of  the 
species  we  can  define  angular  plasma  frequencies  ujvja+’  w^a-’ 

WNm+  anc*  ^Nm  *  Moreover  we  may  employ  a  collisional  frequency 
for  the  electrons  and  a  collision  frequency  v  for  the  ions.  In 
terms  of  these  quantities  we  may  define,  for  an  angular  wave  fre¬ 
quency  w ,  complex  angular  gyro  and  plasma  frequencies  as  follows 
(Booker,  1975): 
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Using  these  quantities,  the  longitudinal,  transverse  and  Hall  suscepti¬ 
bilities  of  a  plasma  at  angular  frequency  u  are  respectively 
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For  vertical  propagation  in  the  ionosphere  at  a  latitude  where  the 
inclination  of  the  Earth's  magnetic  field  to  the  horizontal  is  I,  the 
complex  refractive  indices  of  the  0  and  X  waves  are  given  by 
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The  transition  between  quasi-longitudinal  and  quasi-transverse  pro¬ 
pagation  occurs  where 


2  sin  |  I 
cos  I 


(64) 


At  a  magnetic  pole  of  the  Earth  the  0  and  X  waves  are  the  waves  for 
which  the  electromagnetic  vectors  rotate  respectively  in  the  left 
and  right  handed  senses  round  the  positive  direction  of  the  Earth's 
magnetic  field.  Except  at  the  magnetic  equator  the  0  and  X  waves 
may  also  be  referred  to  as  the  L  and  R  waves  respectively. 

Care  is  necessary  in  evaluating  the  polarization  ratios  Q  used 
in  earlier  sections  because  the  axes  of  coordinates  that  it  has  been 
convenient  to  use  in  this  chapter  differ  from  those  often  used  in  the 
cold  plasma  wave  theory.  The  relation  between  the  sets  of  axes 
involved  is  shown  in  Figure  3.  The  axes  (x,y,z)  are  those  used  in 
this  paper,  with  z  vertically  upwards  and  y  in  the  horizontal  direc¬ 
tion  of  transmission  in  the  Earth-ionosphere  transmission  line. 

Let  us  assume  that  this  direction  points  at  an  angle  A  east  of 
magnetic  north.  The  axes  (N,  W,  V)  in  Figure  3  have  the  respective 
directions  of  magnetic  north,  magnetic  west  and  upward  vertical. 

These  are  the  axes  often  used  by  ionospherists  in  the  cold  plasma 
wave  theory  for  vertical  propagation.  The  V  axis  coincides  with 
the  z  axis  but  the  (N,W)  axes  are  rotated  through  an  angle 
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i  tr  +  A  relative  to  the  (x,y)  axes.  The  polarization  ratio  commonly 
used  in  association  with  the  (N,VJ,V)  axes  is 


P  = 


(65) 


where  E,  and  are  the  westerly  and  northerly  components  of  the 
complex  electric  field  of  the  wave.  The  polarization  ratio  that  it 
has  been  convenient  to  use  in  this  paper  is  (Equations  28  and  29) 
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It  is  easily  shown  that  the  relation  between  the  two  is 


sin  A  +  P  cos  A 
cos  A  -  P  sin  A 


(67) 


However,  it  must  further  be  noticed  that,  if  one  is  to  make  use 
of  the  coefficients  of  susceptibility  listed  in  Equations  (61)  in  a 
susceptibility  tensor  of  the  form 
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the  axes  involved  are  those  marked  (X,Y,Z)  in  Figure  3,  in  which  the  X 
axis  is  westerly  and  the  Z  axis  is  in  the  direction  of  the  Earth's  mag- 
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netic  field.  If  one  performs  the  cold  plasma  wave  calculations  with 
these  axes,  one  uses  the  polarization  ratios  (Booker,  1975) 
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and  in  terms  of  these  expressions  the  polarizat 
from  Figure  3, 
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Sustitution  from  this  equation  into  Equation  (67)  then  gives 

sin  I  sin  A-  Q  cos  A-  Q7  cos  I  sin  A 

q  =  - J : - -  (72) 

sin  I  cos  A+  sin  A-  cos  ^  cos  ^ 

To  obtain  Q0  and  for  use  in  previous  sections  of  this  paper,  we 

2  2  2 
use  for  n  in  Equations  (69)  and  (70)  the  values  n~  and  n^  from 

Equations  (62)  and  then  substitute  into  Equation  (72). 

It  should  be  noted  that  the  limiting  values  of  Q7  in  Equation  (70) 

at  the  equator  (I  =  0)  are 


Q 


2 


i  - 

to 


O  wave 
X  wave 
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This  is  because,  for  Che  0  wave,  the  denominator  of  Q2  also  tends  to 
zero. 
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horizontal  elliptic  polarization  at  ground-level  is  linear  (per¬ 
pendicular  to  the  direction  of  propagation) .  Similar  statements  are 
true  for  the  X  wave  for  equatorial  propagation  in  the  east-west 
direction  (A=  ,  -^it)  . 

We  may  note  that,  if  the  Earth’s  magnetic  field  were  completely 


neglected,  we  would  have 


(78) 


and  Equation  (75)  would  become 


S 


1 


ZRnR 


(79) 


in  accordance  with  Wait  (1962,  p.  290). 

(b)  Quasi-longitudinal  behavior 

Vertical  propagation  of  the  0  and  X  waves  is  quasi-longitudinal 
at  the  level  of  reflexion  in  the  ELF  band  over  a  large  part  of  the 
world.  Under  these  circumstances 


(80) 


(81) 


and  consequently  that 
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These  equations  are  independent  of  A. 

Since  2Q  and  are  normally  not  very  different  numerically, 
we  deduce  from  Equation  (84)  that  the  transmission-line  mode,  under 
quasi-longitudinal  conditions  in  the  ionosphere,  is  nearly  but  not 
quite  linearly  polarized  perpendicular  to  the  direction  of  propagation, 
and  from  Equation  (83)  that  the  degree  of  excitation  of  the  0  and  X 
waves  at  their  respective  levels  of  reflection  is  nearly  but  not  quite 
the  same. 

We  can  also  see  from  Equation  (82)  that,  so  far  as  S  is 
concerned,  the  effects  of  the  0  and  X  waves  are  mathematically  additive. 
However,  the  corresponding  physical  significance  requires  a  rather 
different  statement.  At  the  level  of  ELF  reflection  under  quasi¬ 
longitudinal  conditions  the  0  wave  behaves  largely  evanescently  and 
n^  in  Equation  (82)  is  close  to  a  pure  imaginary,  whereas  the  X  wave 
is  the  whistler  mode  and  n^  in  Equation  (82)  is  close  to  a  real  quantity. 
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This  meant  that  the  0  wave  term  in  Equation  (82)  contributes  mainly 


to  the  real  part  of  S  and  therefore  to  the  departure  of  the  velocity 
of  propagation  in  the  Earth-ionosphere  transmission-line  from  c.  On 
the  other  hand,  the  X  wave  term  in  Equation  (82)  contributes  mainly  to 
the  imaginary  part  of  S  and  therefore  to  the  rate  of  attenuation.  In 
other  words,  the  modification  in  the  velocity  of  transmission  of  ELF 
waves  is  associated  mainly  with  storage  of  energy  by  the  0  wave 
immediately  above  the  level  of  reflection, whereas  attenuation  of  ELF 
waves  in  the  Earth-ionosphere  transmission  line  due  to  leakage  of 
energy  into  the  region  above  the  level  of  reflection  is  associated 
mainly  with  the  X  wave  (the  whistler  wave). 

(c)  Relative  importance  of  the  quasi-longitudinal  and  quasi- 
transverse  approximations 

In  the  ELF  band  there  is  a  drastic  difference  between  the 
practical  importance  of  the  quasi-longitudinal  approximation  and  that 
of  the  quasi-transverse  approximation.  Application  of  Equation  (64) 
shows  that,  for  vertical  transmission  in  the  D  and  E  regions  in  this 
band, propagation  is  quasi-longitudinal  at  practically  all  latitudes. 
The  change  from  quasi-longitudinal  to  quasi-transverse  propagation 
occurs  extremely  close  to  the  magnetic  equator.  Even  at  a  magnetic 
latitude  of  1°,  vertical  propagation  is  quasi-longitudinal  at  all 
heights  less  than  100  km  and  all  frequencies  less  than  300  Hz.  At 
these  frequencies,  therefore,  the  equatorial  zone,  where  there  can  be 
significant  departure  from  quasi-longitudinal  propagation,  is  less 
than  \  /(2:t)  in  width,  and  consequently  has  little  influence  on  world¬ 
wide  propagation  of  ELF  waves  in  the  Earth-ionosphere  transmission 
line,  at  least  by  day. 
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The  almost  complete  dominance  of  the  quasi-longitudinal 


approximation  for  calculation  of  ELF  propagation  in  the  Earth-ionosphere 
transmission  line  has  a  practical  consequence  of  the  greatest  impor¬ 
tance.  Under  quasi-longitudinal  or  quasi-transverse  conditions 
coupling  between  the  0  and  X  waves  is  unimportant.  It  is  only  in  the 
transition  between  quasi-longitudinal  and  quasi-trasnverse  propagation 
that  coupling  is  liable  to  be  important,  and  even  then  not  at  all 
heights.  The  almost  complete  dominance  of  the  quasi-longitudinal 
approximation  implies,  therefore,  that  coupling  between  the  0  and  X 
waves  is  unlikely  to  be  of  much  importance  for  ELF  communications. 

This  constitutes  a  major  simplification  in  the  theory  of  ELF  propa¬ 
gation  in  the  Earth-ionosphere  transmission  line. 
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7. 


The  general  character  of  the  frequency  dependence  of  the  height 
of  reflexion  given  by  Equation  (22)  may  be  illustrated  in  the  following 
way.  Except  very  close  to  the  equator,  vertical  propagation  in  the 
ionosphere  at  ELF  is  quasi-longitudinal.  Moreover,  as  we  move  above 
the  level  where 


v  -  |  x  s in  I  I  (85) 

e  Me  1 

electronic  collisions  contribute  a  relative  minor  correction  to  the 
magnitude  of  the  complex  refractive  index.  Since  the  level  (85) 
occurs  at  about  70  km  except  near  the  equator  and  since  the  level  (22) 
is  frequently  above  the  level  (85) ,  electronic  collisions  usually 
have  only  a  limited  influence  in  determining  the  level  of  reflexion. 

On  the  other  hand,  for  ions  of  angular  gyro  frequency  “Mi,  the 

level 

v  .  -  lx...  sin  I  I  (86) 

i  1  Mi 

occurs  at  about  120  km  or  higher  and  is  appreciably  above  the  level 
of  ELF  reflexion.  Below  the  level  (86)  the  ions  are  prevented  by 
collisions  from  developing  much  motion  and  consequently  from  contribu¬ 
ting  in  a  major  way  to  the  magnitude  of  the  complex  refractive  index. 
Between  the  levels  (85)  and  (86),  whistler-band  behavior  extends  down¬ 
wards  in  frequency  below  the  ionic  gyro-frequency  because  the  ions  are 
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largely  held  at  rest  by  collisions.  Consequently,  over  most  of  the  ELF 
band  and  over  most  of  the  world,  useful  approximations  to  the  refractive 
indices  of  the  0  and  X  waves  are  often  provided  by  the  whistler-band 
formulae 


JU 


X  Me 


sin 


i| 


(87) 


For  an  exponential 
(scale-height  h) ,  Equation 


increase  of  electron  density  with  height 
(87)  has  the  form  of  Equation  (1)  with 


(  0 
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O  wave 
X  wave 
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The  corresponding  values  of  y  in  Figure  1  are 


1  y  (0) 

=  2.  00 

O  wave 

(90) 

f  V  frr ) 

=  2.  57 

X  wave 

(91) 

If  the  quantity 

Me 

(92) 
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evaluated  at  the  reference  height  z  =  a  is  written  u)q  ,  and  if  the 
unity  in  Equation  (87)  is  dropped,  the  height  variations  of  the  refrac¬ 
tive  indices  of  the  0  and  X  waves  become 


2 

n 

O 

2 

n 

X 


so  that 


(93) 


(94) 


Equation  (22)  for  the  heights  of  reflection  zQ  and  zx  of  the  0  and  X 
waves  therefore  gives 


y  (0  )  2tt  h 


Y  (rr)  2rr  h 


(95) 


(96) 


From  these  equations  we  deduce  the  following  approximate  expressions 
for  the  ELF  reflexion  heights: 
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(97) 


1  /  2 

zQ  -  a  +  2h  In  [c  /  fy  (0 )  ( a  uu  )  h  }  ] 

z  =  a  +  2h  In  [c  /  f  y  (tt  )  (  ju  uu  )I/2  hi]  (98) 

X  o  J 


Equations  (97)  and  (98)  illustrate  the  fact  that  heights  of 
reflection  in  the  ELF  band  increase  as  the  frequency  decreases.  This 
behavior  continues  below  the  ionic  gyro- frequency  until  tunnel- 
effect  through  to  the  Alfven  region  above  the  level(86)  results  in 
penetration  of  the  ionosphere  at  micropulsation  frequencies.  The  rise 
in  the  level  of  reflection  as  the  frequency  is  decreased  through  the 
ELF  band,  culminating  in  penetration  of  the  ionosphere  in  the  micro- 
pulsation  band,  is  to  be  compared  and  contrasted  with  the  rise  in 
the  level  of  reflection  as  the  frequency  is  increased  through  the  LF, 

MF  and  HF  bands,  culminating  in  penetration  of  the  ionosphere  in  the 
VHF  band. 

Because  y(0)  <  y(ir)  in  accordance  with  Equations  (90)  and  (91), 
it  follows  from  Equations  (97)  and  (98)  that  Zq  >  z^.  This  illustrates 

the  fact  that  the  level  of  reflection  for  the  0  wave  in  the  ELF  band 
is  normally  a  little  greater  than  that  for  the  X  wave. 

Finally,  let  us  deduce,  for  refractive  index  profiles  based  on 
Equations  (87),  the  propagation  properties  of  the  transmission-line 
mode  in  the  Earth-ionosphere  wave-guide.  From  Equations  (93)-(96),  the 
complex  refractive  indices  for  the  0  and  X  waves  at  their  respective 
levels  of  reflection  may  be  written 
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Substitution  into  Equation  (82)  then  gives 
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V  (0) 
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V  f  T  ) 
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(101) 


For  a  profile  that  is  not  precisely  exponential  we  would  de¬ 
duce  the  scale-height  from  Equation  (23)  and  obtain  slightly 
different  values  hQ  and  hx  at  the  levels  z^  and  z^  respectively. 
Equation  (101)  would  then  become 
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or,  using  Equations  (90)  and  (91)  , 


(102) 
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It  follows  that,  over  a  large  part  of  the  world  under  commonly  en¬ 
countered  ionospheric  conditions,  the  phase  velocity  in  the  Earth- 
ionsophere  transmission  line  is  decreased  below  the  velocity  of  light 
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while  the  attenuation  rate  due  to  leakage  into  the  upper  ionosphere 
in  nepers  per  Xq/(2it)  is  approximately 

h 

0.  64  (104) 

X 


In  decibels  per  megameter  at  a  frequency  f  measured  in  hertz,  the 
attenuation  rate  is  approximately 


h 

0.  1 17  f  -X 
z 

X 


(105) 
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Chapter  3.  Theory  of  ELF  Propagation  in  the  Earth- 
ionosphere  Transmission  Line 


1.  Introduction 

Numerous  attempts  have  been  made  to  devise  simple  methods 
for  handling  ELF  propagation  of  the  transmission-line  mode  of  propa¬ 
gation  between  the  conducting  Earth  and  the  ionosphere.  Many  of  these 
are  summarized  in  books  by  Budden  (1961a, b).  Wait  (1962)  and  Galejs 
(1972),  and  also  in  a  special  issue  of  the  IEEE  Transactions  on 
Communications  edited  by  Wait  (1974).  Further  efforts  in  this  direc¬ 
tion  have  been  made  by  Booker  and  Lefeuvre  (1977)  as  described  in 
Chapter  2,  and  by  Greifinger  and  Greifinger  (1978,1979).  These  pairs 
of  authors  will  be  referred  to  as  BL  and  GG. 

Both  BL  and  GG  find  that  there  is  an  identifiable  reflecting 
stratum  for  the  return  of  ELF  waves  from  the  ionosphere.  Moreover, 
in  spite  of  different  approaches,  they  come  to  the  same  conclusion 
concerning  the  location  of  the  stratum.  The  thickness  is  about  one 
scale-height  h,  defined  in  terms  of  the  profile  of  n^-l,  where  n  is 
the  ionospheric  refractive  index.  The  center  of  the  reflecting 
stratum  occurs  where  1/(2tt)  times  the  local  wavelength  \  is  equal  to 
Y  times  the  local  scale-height  h,  where  y  is  a  number  between  2  and  3. 
The  BL  treatment  of  the  value  of  y  is  slightly  different  from 
that  of  GG.  Preference  will  be  given  to  the  BL  treatment. 

Above  the  reflecting  stratum,  the  ionosphere  is  a  good  approxi¬ 
mation  to  a  slowly  varying  medium.  It  is  a  slowly  varying  medium  in 
which  the  ordinary  (0)  and  extraordinary  (X)  waves  can  nearly 
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always  be  assumed  Co  propagate  independently.  Moreover,  because  of  the 
high  numerical  value  of  the  refractive  index,  phase  propagation  above 
Che  reflecting  stratum  (but  not,  in  general,  group  propagation)  is 
almost  exactly  vertical  regardless  of  the  directions  of  propagation 
of  the  incident  and  reflected  waves  below  the  ionosphere.  Major 
simplifications  arise  from  these  considerations. 

Neither  BL  nor  GG  draw  explicit  attention  to  another  important 
feature  of  the  reflecting  stratum.  When  a  downcoming  magneto-ionic 
component  (0  wave  or  X  wave)  leaves  the  ionosphere,  a  location  is 
reached  below  which  no  further  important  change  in  elliptical  polari¬ 
zation  takes  place.  This  occurs  because  a  magnetoplasma  does  not  react 
in  a  doubly  refracting  manner  unless  the  refractive  indices  of  0  and 
X  waves  differ  by  more  than  about  A/(27rh)  where  A  is  the  local  wave¬ 
length  and  h  is  the  local  scale-height  (Booker  1936,  Budden  1952,  Hayes 
1971).  In  the  HF  band  the  level  thus  determined  can  usually  be  re¬ 
garded  as  the  "bottom'1  of  the  ionosphere.  In  the  ELF  band,  however, 
because  of  the  much  bigger  free-space  wavelength  A  ,  the  level  below 
whic?  the  ionosphere  fails  to  respond  as  a  doubly  refracting  medium 
is  approximately  the  same  as  the  level  of  reflection.  Consequently 
the  reflecting  stratum  as  identified  by  BL  and  GG  has  the  property 
that,  above  it,  the  ionosphere  is  a  slowly-varying  non-isotropic 
medium  whereas,  below  it,  the  ionosphere  responds  primarily  as  an 
isotropic  medium. 

However,  BL  overlooked  an  important  feature  taken  into  account 
by  GG.  BL  took  the  isotropic  medium  below  the  reflecting  stratum  as 
free  space,  allowing  for  collisional  absorption  in  this  region  by  a 
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method  similar  to  that  employed  for  dielectric  loss  in  an  engineering 
transmission  line.  GG  emphasized  that  the  isotropic  medium  below  iue 
level  of  reflection  is  strictly  equivalent  to  free  space  only  below  the 
level  where  the  refractive  index  begins  to  approach  unity.  Between  this 
level  and  the  level  of  reflection,  propagation  is  essentially  vertical 
almost  regardless  of  the  directions  of  incidence  and  reflection  below 
the  bottom  of  the  ionosphere.  At  ELF  the  transition  from  nearly 
horizontal  propagation  to  nearly  vertical  propagation  takes  place,  not 
in  the  reflecting  stratum,  but  in  a  refracting  stratum  at  the  bottom 
of  the  ionosphere.  What  GG  show  is  that  the  vertical  standing-wave 
pattern  for  the  horizontal  electric  field  below  the  level  of  reflection 
is  substantially  what  it  is  in  free  space  subject  to  the  following 
remarkable  proviso.  Propagation  of  phase  is  in  the  nearly  horizontal 
directions  of  incidence  and  reflection  below  the  bottom  of  the  ionosphere 
but,  between  the  bottom  of  the  ionosphere  and  the  level  of  reflection, 
change  of  phase  is  to  be  calculated  to  a  first  approximation  on  the 
basis  of  vertical  propagation  in  free  space. 

This  unusual  result  is  an  indispensable  key  to  understanding  ELF 
propagation.  The  way  in  which  the  various  components  of  the  electro¬ 
magnetic  field  behave  is  described  by  GG  using  wave  theory. 

The  result  seems  unusual  because,  between  the  refracting  stratum  at  the 
bottom  of  the  ionosphere  and  the  reflecting  stratum  higher  up,  we  are 
dealing  with  a  rapidly  varying  stratifed  medium  and  there  is  a  lack  of 
familiarity  with  solutions  of  the  wave  equation  in  these  circumstances. 

In  particular,  it  is  essential  to  realize  that  there  is  remarkably 
little  change  of  phase  in  passing  through  such  a  rapidly  varying 
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medium  in  comparison  with  what  would  be  mistakenly  calculated  using 


the  theory  for  a  slowly  varying  medium. 

In  this  chapter  the  ideas  of  BL  and  GG  are  combined  into  a  single 
approach.  In  taking  the  Earth's  magnetic  field  into  account  there  is  no 
restriction  to  quasi-longitudinal  propagation.  In  Chapter  4  the  theory 
is  applied  to  a  simple  worldwide  model  of  the  ionosphere. 
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!.  The  equivalent  complex  height  of  reflection 

We  neglect  the  curvature  of  the  Earth  and  suppose  that  the 
horizontal  direction  of  propagation  in  the  Earth-ionosphere  trans¬ 
mission  line  is  at  an  angle  A  east  of  north.  We  take  the  magnitude 
and  direction  of  the  Earth's  magnetic  field  to  be  specified  functions 
of  the  height  z  above  the  Earth's  surface.  Likewise  we  take  the 
electron  density,  the  ion  densities  and  the  collisional  frequencies 
to  be  specified  functions  of  z. 

Although  propagation  of  phase  is  nearly  horizontal  below  the 

ionosphere,  we  take  it  to  be  nearly  vertical  above  the  bottom  of  the 

ionosphere.  For  vertical  propagation,  calculations  can  be  made  of 

the  complex  refractive  indices  n^  and  n^  for  the  0  and  X  waves,  as 

well  as  their  complex  elliptical  polarizations  Oq  and  Q^.  The 

necessary  formulae  are  summarized  in  Section  5  of  Chapter  2,  and  a  more 

complete  account  is  given  by  Booker  (1975).  It  should  be  noted  that 

the  definitions  of  the  0  and  X  waves  used  by  BL  may  not  be  the  most 

convenient.  There  exist  a  height  zr  and  a  magnetic  latitude  A  where 

L  C 

there  is  a  branch  point  at  which  n^  •  n^.  The  branch  point  is 
located  in  the  E  region  close  to  the  magnetic  equator.  BL  maintained 
continuity  of  wave  nomenclature  at  latitude  for  z  >  z^,  whereas  it 
might  be  more  convenient  to  do  so  for  z  <  z^.  Nevertheless  we  shall 
preserve  the  BL  nomenclature. 

At  each  height  z  in  the  ionosphere  we  define  the  scale  heights 
1)q  and  for  vertically  propagating  0  and  X  waves  by  means  of  the 
equations  (Equation  23  of  Chapter  2) 
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1  _d_  ,  2  "1  _  l  JL  fn2  n 

2  .  dz  (n0  _1)  *  hX  2  dz  (nX  _1) 

nrt  -1  Ov 


The  levels  of  reflection  z^  and  z^  for  the  0  and  X  waves  are  then 
such  that  (Equation  22  of  Chapter  2) 


x  2rThx 


where  the  numbers  y.  and  y^  are  read  from  Figure  1  of  Chapter  2. 

The  procedure  for  determining  y^  and  y^  is  one  of  successive  approxi¬ 
mation.  First  assume  that  yQ  =  2.00,  Yx  =  2.57  and  use  Equations  (2) 
to  determine  a  first  approximation  to  zQ  and  z^.  For  these  heights 
evaluate  (Equation  24  of  Chapter  2) 


=  T  +  2  arg  n. 


=  tt  +  2  arg  n x 


Use  these  values  of  and  in  Figure  1  of  Chapter  2  to  determine 
improved  values  of  y^  and  y^.  Use  these  values  in  Equations  (2)  to 
obtain  improved  values  for  Zq  and  z^. 

Knowing  the  levels  of  reflection  z^  and  z^,  and  the  correspond¬ 
ing  values  n_  and  n  of  the  refractive  indices,  we  are  now  in  a 
position  to  calculate  equivalent  complex  heights  of  reflection  for 
the  0  and  X  waves.  These  are  defined  by  the  equations 


^RO  "  Z0  “  J  27in  ’  ZRX  =  ZX  “  J  2^ 
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In  Section  4  these  equivalent  complex  heights  of  reflection  will  turn 
out  to  play  a  key  role  in  calculating  the  phase  velocity  and  attenua¬ 
tion  rate  in  Earth- ionosphere  transmission  line. 

It  should  be  noted  that,  because  n^  and  n^  are  complex,  the 

real  parts  of  z '  and  ~  differ  from  the  actual  heights  of  reflec- 
RU  Ra 

tion  z .  and  z  .  Except  near  the  equator,  it  is  quite  common  for  n 

U  A  A 

to  be  almost  real  and  for  n^  to  be  almost  purely  imaginary.  The  real 

part  of  z  in  the  second  of  Equations  (4)  is  then  almost  equal  to  z  , 

and  there  is  an  imaginary  part  linked  primarily  to  leakage  of  energy 

into  the  upper  ionosphere  via  the  whistler  mode.  But  the  transmitted 

0  wave  behaves  evanescently .  Consequently  z^q  in  the  first  of 

Equations  (4)  is  nearly  real  but  larger  than  z^;  this  is  linked  to  a 

phase  change  on  reflection  that  affects  the  velocity  of  propagation  in 

the  Earth-ionosphere  transmission  line.  In  general,  there  are 

negative  imaginary  parts  to  z' _  and  z’  ,  and  also  real  parts  that  differ 

KU  Ra 

from  Zq  and  z^  respectively.  It  is  frequently  true  that  the 

difference  between  zn  and  z  is  small,  but  it  is  seldom  feasible  to 

neglect  the  difference  between  z"  _  and  Z*  . 

RU  RX 
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3.  The  equivalent  complex  height  of  the  bottom  of  the  ionosphere 

The  GG  method  for  taking  into  account  the  effect  of  ioniza¬ 
tion  below  the  level  z  of  reflection  depends  on  the  assumption  that 

K 

the  horizontal  complex  magnetic  field  H(z)  of  the  wave  in  the  Earth- 

ionosphere  transmission  line  is  almost  independent  of  height  between 

z  =  0  and  z  =  z  .  For  an  angle  of  incidence  whose  sine  is  S (complex 
R 

for  a  wave-guide  mode),  the  wave  equation  for  H(z)  is 


(5) 


Let  E(z)  be  the  corresponding  horizontal  component  of  the  complex 
electric  field.  E(z)  is  derived  from  H(z)  by  means  of  the  equation 


E(z) 


o 

"J  2 IT 


n2(zr) 


dH 

dz 


(6) 


where  is  the  impedance  of  free  space.  By  integrating  Equation  (5) 
with  respect  to  z  and  substituting  into  Equation  (6)  it  follows  that, 
if  the  horizontal  electric  field  vanishes  at  ground-level,  then  at  the 
level  of  reflection  we  have 


E(ZR} 


j  C, 


2tt 

> 

o 


s2  ) 

—7 -  H(z)dz 

n  (z)  ' 


(7) 


The  assumption  is  now  made  that,  in  this  integral,  H(z)  may  be  taken 
as  constant  across  the  Earth-ionosphere  transmission  line.  It  follows 
that 
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E(zr) 


^oH 


2t 


n2(z)  ) 


dz 


(8) 


and  therefore  that 

E(V  -  j  «0«’/y(--R-s2-?B>H 

where 


(9) 


(10) 


The  integrand  in  Equation  (10)  is  unity  below  the  ionosphere 

but,  in  the  ELF  band,  it  decreases  rapidly  as  the  ionosphere  is  entered. 

Consequently  ‘z  is  a  measure  of  the  height  of  the  bottom  of  the  ionos- 
B 

phere.  However,  because  n(z)  is  complex,  it  follows  that  t  is  complex. 
As  shown  by  GG,  the  positive  imaginary  part  of  zl  plays  an  important 

D 

role  in  connection  with  attenuation  in  the  Earth-ionosphere  transmission 

line  caused  by  ionization  below  the  level  of  reflection.  The 

quantity  (s  is  the  equivalent  complex  height  of  the  bottom  of  the 
B 

ionosphere. 

Equation  (9)  arises  from  the  vertical  standing-wave  pattern 
existing  between  the  Earth's  surface  and  the  level  of  reflection  in 
the  ionosphere.  The  quantity  S  is  the  sine  of  the  complex  angle  of 
incidence  below  the  ionosphere.  When  the  boundary  conditions  at 
the  level  of  reflection  are  satisfied,  S  is  the  eigenvalue  of  the 
mode  problem.  Its  real  and  imaginary  parts  give  the  propagation 
properties  of  the  Earth-ionosphere  transmission  line. 
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2  2 

In  Equation  (9)  let  us  replace  S  by  1-C  ,  so  that  C  is  the 
cosine  of  the  complex  angle  of  incidence  upon  the  ionosphere.  We  then 
obtain 


(  2ttz  -  2tt(z  -  z  ) 

e(zr>  -  + — x — 

l  0  o 


(11) 


The  corresponding  result  in  the  BL  treatment  of  the  region  below  the 
# 

level  of  reflection  is  (Section  3,  Chapter  2) 


e(zr> 


J  5, 


(12) 


Equation  (12)  arises  from  the  vertical  standing-wave  pattern  of 
horizontal  electric  field  existing  between  the  Earth's  surface  and  the 
level  of  reflection  on  the  assumption  that  we  have  waves  propagating 
in  free  space  at  an  angle  of  incidence  and  an  angle  of  reflection 
whose  cosine  is  C.  Equation  (11)  shows  that  this  is  only  appropriate 
below  the  bottom  of  the  ionosphere;  between  the  bottom  of  the  iono¬ 
sphere  and  the  level  of  reflection,  the  vertical  standing-wave  pattern 
for  horizontal  electric  field  should  be  taken  to  be  that  for  vertical 
propagation  in  free  space.  What  the  work  of  GG  shows  is  that  the  BL 
treatment  of  the  standing  wave  behavior  across  the  Earth-ionosphere 
transmission  line  must  be  modified  by  replacing  the  right  hand  side  of 
Equation  (12)  by  the  right  hand  side  of  Equation  (11).  This  is  easy 
to  do  and  the  results  are  given  in  the  following  section. 

It  may  be  noted  that  two  definitions  are  possible  for  the 
equivalent  complex  height  of  the  bottom  of  the  ionosphere  depending 
on  whether  the  complex  refractive  index  in  Equation  (10)  is  taken  to 
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be  that  for  the  0  wave  or  that  for  the  X  wave.  The  two  definitions 


are : 


There  is  not  much  difference  between  the  two  because  n^  and  n^  only 
differ  where  they  are  numerically  large.  However,  for  north-south 
transmission  at  the  equator,  only  the  0  wave  is  excited  and  "z ’  is  to 
be  preferred  while,  for  east-west  transmission  at  the  equator,  only  the 
X  wave  is  excited  and  z'  is  to  be  preferred.  On  the  other  hand,  at  a 
magnetic  pole,  transmission  in  any  horizontal  direction  excites  the  0 
and  X  waves  above  the  level  of  reflection  almost  equally  and  an 
average  of  z^  and  ‘Z’  is  to  be  preferred.  It  is  appropriate  to  weight 
and  ?BX  according  to  the  relative  degree  of  excitation  of  the  0 
and  X  waves.  This  may  be  conveniently  done  in  the  following  way. 

As  described  in  Chapter  2,  the  mode  in  the  Earth-ionosphere 
transmission  line  consists  of  a  quadruplet  of  crossing  waves,  one 
pair  having  the  elliptical  polarization  of  the  0  wave  just  above  its 
level  of  reflection,  and  the  other  pair  having  the  elliptical 
polarization  of  the  X  wave  just  above  its  level  of  reflection.  For 
the  pair  of  crossing  waves  with  0  polarization  we  shall  write  the 
factor  in  curly  brackets  on  the  right  hand  side  of  Equation  (11)  as 


+ 


(14) 
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where  z"  is  given  by  the  first  of  Equations  (13) .  On  the  other  hand 
BU 

for  the  pair  of  crossing  waves  with  X  polarization  we  shall  write 


2lTZBX  „2  ,  2TT(zX_ZBx2 

\  L  r 


(15) 


X. 


where  z'  v  is  given  by  the  second  of  Equations  (13) 
BX 
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4 .  Propagation  velocity  and  attenuation  rate 


We  are  now  in  a  position  to  adapt  the  BL  theory  of  propaga¬ 
tion  in  the  Earth-ionosphere  transmission  line  to  incorporate  the  GG 

theory  of  the  effect  of  ionization  below  the  level  of  reflection.  It 

2 

is  a  matter  of  replacing  the  respective  quantities  (2itZq/Xo)C  and 
2 

(2ttzx/Ao)C  in  Equations  (38)  -  (41)  of  BL  by  expressions  (14)  and 
(15)  in  this  Chapter,  and  then  ascertaining  the  consequent  changes  in 
Equations  (46)  -  (48)  of  Chapter  2.  The  upshot  is  that  Equation  (51) 
of  Chapter  2  may  be  written 


S 


,-l. 


Q0  ^B0^R(P  _QX  ^BX^RX^ 


(16) 


where  z"RQ  and  “z^  are  given  by  Equations  (4),  2"^  and  2^  are  given  by 
Equations  (13)  ,  and  Qq  and  are  the  complex  elliptical  polarizations 
of  the  0  and  X  waves  at  their  respective  levels  of  reflection.  When 
the  complex  eigenvalue  S  has  been  calculated  from  Equation  (16)  ,  the 
ratio  of  the  velocity  c  of  light  to  the  phase  velocity  v  in  the  Earth- 
ionosphere  transmission  line  is  given  by 

c/v  -  R (S)  (17) 

and  the  rate  of  attenuation  in  nepers  per  unit  distance  by 

a  =  -(2tt/Xq)  I  (S) .  (13) 

where  R  and  I  denote  real  and  imaginary  parts. 

Over  a  substantial  part  of  the  Earth,  vertical  propagation  in 
the  ionosphere  is  quasi-longitudinal  at  ELF  and  we  have  (Equations  81 
of  Chapter  2) 
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Equation  (16)  then  becomes 

S  -  lit  C^q^q>  +  ^BX^RX} }  1 
If  we  further  assume  that 

^  ,  -v  .  ^ 

ZB0  "  ZBX  ?  ZB 

Equation  (20)  becomes 


S  - 


(19) 


(20) 


(21) 


(22) 


If  in  addition  we  assume  that 

^  ,  /V 

"RO  •  ZRX  •  ZR 
Equation  (22)  simplifies  to 

s  ■ 

This  corresponds  to  a  result  derived  by  GG.  However,  while  the  approxi¬ 
mations  in  Equations  (19)  and  (21)  are  often  quite  good,  the  approxi¬ 
mation  in  Equation  (23)  is  inadvisable.  It  is  true  that  the  actual 
heights  of  reflection  of  the  0  and  X  waves  are  often  nearly  the  same, 
so  that 


(23) 


(24) 
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But  the  corrections  required  in  calculating  z^q  and  in  accordance 
with  Equations  (4)  are  significantly  different,  because  the  X  wave 
is  the  whistler  mode  (n^  almost  real)  whereas  the  0  wave  is  quite 
evanescent  (n^  almost  purely  imaginary) . 

Simplification  of  Equation  (16)  also  occurs  at  the  magnetic 
equator.  For  propagation  in  the  Earth- ionosphere  transmission  line 
at  an  angle  A  east  of  north  we  have  at  the  magnetic  equator  (Equations 
74  of  Chapter  2) 

Qq  -  tan  A,  Qx  -  -  cot  A  (26) 

so  that  Equation  (16)  becomes 

S"  C^B0^R0)cOs2a  +  (27) 

For  north-south  transmission  (A  =  0,  tt)  we  have 

s  ■  (:8> 

and  for  east-west  transmission  (A  =  t  4  tt)  we  have 

*  •  (W*  (29> 

In  principle,  propagation  in  the  Earth-ionosphere  transmission  line 
at  the  magnetic  equator  is  a  function  of  the  direction  A  of  horizontal 
transmission,  with  only  the  0  wave  excited  above  the  level  of 
reflection  for  north-south  transmission  and  only  the  X  wave  for  east- 
west  transmission.  In  practice,  the  numerical  values  of  the  ex¬ 
pressions  on  the  right  hand  sides  of  Equations  (28)  and  (29)  are 
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nearly  equal  at  the  equator,  and  no  significant  dependence  on  the 
angle  of  transmission  A  materializes.  This  is  because  vertical  propa¬ 
gation  of  both  the  0  and  X  waves  immediately  above  the  level  of 
reflection  at  the  magnetic  equator  is  largely  metallic,  so  that  n^  and 
n  in  Equations  (4)  are  nearly  equal  both  in  absolute  value  and 

A 

argument,  the  latter  being  close  to  -f/4. 

It  may  be  noticed  that,  in  Equations  (16),  (20)  and  (27),  the 
imaginary  parts  of  the  complex  heights  are  usually  small  compared 
with  the  real  parts,  and  that  useful  binomial  expansions  may  be  per¬ 
formed  on  this  basis. 

The  approximations  discussed  in  this  section  are  often  convenient 
for  visualizing  the  physics  of  the  propagation  phenomenon.  These 
approximations  are,  however,  unnecessary  in  numerical  calculations. 
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5 .  Simultaneous  reflection  from  the  D  and  E  regions 

So  far,  we  have  assumed  that  only  a  single  reflecting  stratum 

exists  for  each  magneto-ionic  component  (0  wave  or  X  wave) .  This 

situation  is  illustrated  in  Figure  1.  For  the  selected  magneto-ionic 

component,  the  relevant  Equation  (2)  has  the  solution  z  =  z  .  At 

this  level  the  complex  refractive  index  for  the  magneto-ionic  component 

is  n  .  Reflection  is  calculated  as  though  from  a  discontinuity  of 
K 

refractive  index  n_-l  at  this  height.  Below  this  level  the  variation 
of  complex  phase  with  height  is  calculated  by  the  theory  of  GG,  and 
this  involves  a  refracting  stratum  near  the  bottom  of  the  ionosphere 
where  nearly  vertical  propagation  is  converted  to  nearly  horizontal 
propagation.  The  calculation  based  on  an  individual  magneto-ionic 
component  does  not  satisfy  the  boundary  conditions  at  the  surface 
of  the  Earth  because  the  wave  is  ellipticallv  polarized,  but  the  combi¬ 
nation  of  the  0  and  X  wave  calculations  does  satisfy  the  boundary 
conditions.  Below  the  reflecting  stratum  the  medium  behaves  sub¬ 
stantially  isotropically;  it  is  the  combined  field  derived  from  the 
0  and  X  wave  calculations  that  is  then  significant. 

However,  it  is  quite  easy  for  more  than  one  reflecting  stratum 
to  exist  below  the  level  of  maximum  ionization  density  in  the  E 
region.  This  occurs  when  the  profile  of  ionization  density  has  a  high 
gradient  on  the  under  side  of  the  D  region  and  a  high  gradient  on  the 
under  side  of  the  E  region,  with  relatively  low  gradients  between  and 
above.  As  we  shall  see  in  detail  in  the  next  chapter,  the  behavior 
in  the  ELF  band  is  then  as  follows. 

The  primary  reflecting  stratum  at  the  Schumann  resonant 
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Chapter  3 

Figure  1.  Illustrating  reflection  of  the  0  or  X  wave  from  the 
ionosphere  at  ELF  when  the  relevant  Equation  (2)  has 
only  one  solution. 


frequency  ( —  7.5  Hz)  is  located  in  the  high  gradient  on  the  under  side 
of  the  E  region.  But,  as  the  frequency  increases  through  about  30  Hz, 
partial  reflection  from  the  D  region  comes  into  play.  Below  about 
100  Hz,  the  D  region  is  an  optically  thin  stratum  that  is  quite 
transparent,  so  that  reflection  is  still  largely  from  the  £  region. 

But  with  sufficient  increase  of  frequency,  the  high  gradient  on  the 
under  side  of  the  D  region  takes  over  as  the  primary  seat  of  reflection. 

This  behavior  manifests  itself  in  the  process  of  solving 
Equations  (2)  to  obtain  the  heights  of  reflection  Zq  and  z^.  For  the 
two  waves,  let  us  introduce  local  wavelengths  Aq  and  A^  defined  in 
terms  of  the  wavelength  X  in  free  space  by 

-  VKt’  Xx”  (30) 

Equations  (2)  may  then  be  written 


By  plotting  the  quantities  on  the  left-hand  sides  of  these  equations 
as  functions  of  height,  it  is  easy  to  locate  the  heights  of  reflection 
in  the  manner  shown  in  Figure  2.  Using  a  model  of  the  ionosphere  of 
the  type  employed  in  the  next  chapter.  Figure  2  is  drawn  for  day¬ 
time  conditions  in  a  middle  lattitude.  The  two  functions  of  height 
on  the  left-hand  sides  of  Equations  (31)  are  so  nearly  identical  that 
duplicate  curves  cannot  easily  be  shown  for  the  0  and  X  waves  in 
Figure  2.  The  four  curves  refer  to  four  frequencies.  The  first 
approximations  to  and  z^  are  given  by  the  intersections  of  these 
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2-  Illustrating  the  determination  of  the  day-time  heights 
of  reflection  at  four  frequencies. 
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Figure  3.  Illustrating  multiple  solutions  of  Equations  (2)  under 
day-time  conditions. 
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curves  with  the  ordinates  at  2.00  (0  wave)  and  2.57  (X  wave).  After 


correction  by  successive  approximation  as  described  in  Section  2, 
the  values  of  zn  and  z  behave  as  showr  in  Figure  3. 

For  each  wave,  there  is  only  a  single  solution  of  Equations  (2) 
when  the  frequency  is  sufficiently  low  or  sufficiently  high.  But 
there  is  an  intermediate  band  of  frequencies  where  three  solutions 
exist.  In  this  intermediate  band  the  ionosphere  is  rendered  isotropic 
both  by  the  high  gradient  on  the  under  side  of  the  D  region  and  by 
the  high  gradient  on  the  underside  of  the  E  region.  But  the  ionosphere 
functions  as  a  slowly  varying  non-isotropic  medium  near  the  maximum 
of  ionization  density  in  the  D  region  and  also  near  the  maximum  of 
ionization  density  in  the  E  region.  We  shall  idealize  this  situation 
by  applying  the  phase-integral  treatment  in  the  region  of  Figure  3 
marked  as  slowly  varying  and  non-isotropic,  and  a  GG  type  solution  in 
the  region  marked  as  rapidly  varying  and  isotropic;  the  latter  becomes 
the  solution  appropriate  to  a  homogeneous  isotropic  medium  below  the 
refracting  stratum  at  the  bottom  of  the  ionosphere. 
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6.  Theory  for  three  reflection  heights 

Consider  a  situation  in  which,  for  a  particular  magne' o- 
ionic  component  (0  wave  or  X  wave) ,  there  are  three  heights  of  re¬ 
flection  z^,  z 2  and  z  ( z ^  <_  z^  <_  z^)  as  illustrated  in  Figure  3.  Let 
n^,  n^  and  n^  be  the  complex  refractive  indices  for  this  wave  at  these 
heights.  Let  Q^,  and  be  the  corresponding  complex  elliptical 
polarizations.  While  n^,  n^  and  n^  are  always  markedly  different,  we 
have  found  in  our  calculations  that  the  elliptical  polarizations  at 
the  three  levels  are  nearly  always  closely  the  same.  We  shall  there¬ 
fore  assume  that 

Q1  "  Q2  “  Q3  "  Q  (32) 

but  that 

I4  n2  ^  n3  (33) 

We  use  the  axes  of  coordinates  (x,  y,  z)  shown  in  Figure  3 
of  Chapter  2,  taking  the  positive  y  axis  to  be  the  horizontal 
direction  of  transmission  in  the  Earth-ionosphere  transmission  line 
and  the  positive  z  axis  to  be  vertically  upwards.  For  a  vertically 
propagating  magnetic-ionic  component,  the  complex  phase-change  between 
height  z^  and  height  z ^  is,  according  to  the  phase-integral  theory, 

Z2 

n  dz  (34) 

1 

and  that  between  height  z^  and  height  z^  is,  according  to  the  GG 
theory, 

*23  “  (2ll/Xo)(z3"Z2)  (35) 


12 


(2tt/X  ) 
0 
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to  those  at  height  z^  is  (cf  Equation  9) 


K0(2»/X0)(2l-S%> 


(40) 


where 


7  -  P-O- 

B  1  n2(z) 


(41) 


By  multiplying  the  matrix  (38)  by  the  matrix  (40)  we  see  that  the  matrix 
relating  the  values  of  E  and  H  at  the  Earth's  surface  to  those  at 

XV 

height  is 


a  -  c(27r /A  ) z 
o  1 

LjCo_1{a(2TrAo)(z1-S2IB)+c} 


j  ;  {b  +  d(2irA „)*.}  . 

o  1  1  (42) 

d-b(2rrA0)(z1-S^SB)/ 


We  also  require  corresponding  matrices  for  E^  and  H  .  If 

these  are  written  as  matrices  for  E  and  -H  ,  the  matrix  (38)  is 

y  x 

unchanged  and  the  matrix  (40)  is  reflected  in  the  leading  diagonal.  It 

follows  that  the  matrix  relating  the  values  of  E  and  H  at  the  Earth's 

y  x 

surface  to  those  at  height  is 


a  -  c(2ttAq)  (z1  -  S2^) 
-j  C  "1(a(2-rr/Xo)z1+  c} 


2~ 


-j  ;  {b  +  d(2ir/\o)(z1  -S  z0) 


(43) 


d  -b(27T/X  )z. 

O  1 


With  the  aid  of  the  matrices  (42)  and  (43),  we  are  in  a 

position  to  express  the  horizontal  components  E  ,  H  ,  E  ,  H  of  the 

x  y  y  x 
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complex  electromagnetic  field  at  the  Earth’s  surface  in  terms  of  the 
horizontal  components  E^,  Hy3,  Ey3’  Hx3  at  ttie  height  Zy  Moreover, 
at  the  height  z^,  we  have  (cf  Equations  26-29  of  Chapter  2) 


E  ,  E  , 

-Zl  =  =  _L  r 

H  .  H  «  n,  Sd 

x3  y3  3 


E  ,  H  , 

-Si  .  .  Q 

E  .  U  ,  4 

y3  x3 


These  equations  permit  us  to  express  E^,  H^3  and  E^3  in  terms  of 

H  If  we  write 

x3 


f  ■  (a/n.,)  +  j  b. 


g  *  d  +  j(c/n3) 


we  obtain 


>o  Ex  -  U  +  jg(2TT/Xo)Zl}QHx3 

Hy  =  {j  f  (2rr/Ao)  <Z;L  -  s\)  +  g}Q  Hx3 
C^y  =  -tf  +  j8(27r/Ao)(Zl-S27B)}Hx3 
Hx  “  (j  f(2TT/Ao)Zl  +  g}Hx3 


In  the  second  and  fourth  of  Equations  (47)  the  g  terms  dominate, 
and  these  equations  may  be  simplified  to 


V  8  « »,3  f 
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But  in  the  first  and  third  of  Equations  (47)  it  is  necessary  to  avoid 

approximations.  However,  an  equivalent  complex  height  of  reflection 

z'_  may  be  introduced  such  that 
K 

f  +  j  g(2TT/Xo)z1  =  j  g(27T/XQ)7R  (49) 

The  first  and  third  of  Equations  (47)  then  become 


=  jg(2lt/V\QHx3  | 

-  -j  g(2*A0)(\-s\)Hx3  j 


With  appropriate  affixes  0  and  X,  Equations  (48)  and  (50) 
apply  to  both  the  0  wave  and  the  X  wave.  Taking  both  waves  into 
account,  the  complete  expressions  for  the  horizontal  components  (H^, 

Hy,  E^,  Ey)  of  the  complex  electromagnetic  field  at  the  Earth's  surface 
are  given  by 


“x  -  «0  »X3  +  %  "*3  j 
Hjr  "  g0  Q0  “x3+  8XQxHx3  I 


(51) 


and 


-•J  —  —  E 
3  2*  C  * 

X  1 

1  —  —  E 
3  2tt  C  y 

o 


g0  ZR0  Q0  Hx3 


+  gX  ZRX  QX  Hx3 


(52) 


g0^R0  "  S%0)Hx3  +  gX^RX  ‘  s2zBX)Hx3 
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We  are  now  in  a  position  to  satisfy  the  boundary  condi¬ 
tions  at  the  surface  of  the  Earth.  Leakage  of  energy  into  the  Earth 
is  easily  allowed  for  separately  so  that,  in  studying  the  effect  of 
the  ionosphere,  we  may  take  the  boundary  conditions  at  the  surface 
of  the  Earth  to  be 

E  -  E  -  0  (53) 

x  y 

Substitution  from  Equations  (53)  into  Equations  (52)  gives 


Q0  ZR0  80  Hx3  +  QX  2RX  8X  Hx3*  ° 

“S  zB0)g0  Hx3+  ^RX  “S  ZBX)sXHx3 


(54) 


ox 

and  elimination  of  and  then  recovers  Equation  (16) .  The 

relative  degree  of  excitation  of  the  0  and  X  waves  above  the  height 
Zy  is  given  by  the  first  of  Equations  (54)  while  the  horizontal  magnetic 
field  at  the  Earth's  surface  is  given  by  Equations  (51). 

We  see  that,  when  Equations  (2)  for  the  heights  of  reflection 
of  the  0  and  X  waves _each  have  a  triplet  of  solutions  instead  of  a 
single  solution,  the  procedure  for  calculating  the  phase  velocity  and 
rate  of  attenuation  in  the  Earth-ionsosphere  transmission  line  is 
affected  in  only  one  respect.  The  equivalent  complex  heights  of 
reflection  for  the  0  and  X  waves  must  be  calculated,  not  from 
Equations  (4),  but  from  Equation  (49).  This  equation  gives 

v  *!-*»!  (55) 
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T 


or,  using  Equations  (46)  and  (39), 

~  \  n2  B  -  1A 

ZR  “  21  2tt  n1n;}  D  +  jC 

where 

A  -  cos  4>12  cos  <$23  -  (l/n2)  sin  $12  sin  ^  ' 

B  *  n^  cos  <J>32  sin  423  +  (n^/n^  sin  $^2  cos  $23 

C  -  (n2/n3)  sin  0^  cos  $23  +  (l/n3)  cos  <>12  sin  <J>23 
D  -  -  n2  sin  <J>12  sin  <J>23  +  cos  ^  cos  $23  / 


(56) 


(57) 


The  complex  phase-change  $23  in  these  equations  is  always 
small,  and  so  are  the  terms  involving  l/n2  and  l/n3>  Hence  Equations 
(57)  may  usually  be  simplified,  if  desired,  to 

A  ■  cos  <t>12  \ 

B  -  n3  $23  cos  <f^2 +'(n3/n2)  sin  $12  / 

i  (58) 

C  •  (n2/n3)  sin  1>12 
D  ■  -  n2  <t23  sin  $^2  +  cos  $^2  , 


At  the  low-frequency  end  of  the  ELF  band  where  reflection  from  the  D 
region  is  comparatively  unimportant,  the  complex  phase-change  $^2  is 
also  small  and  Equations  (58)  become  approximate!-.- 


A  -  1 


B 

C 


n3  $23  +  (n3/n2)  *12 
(n2/n3)  *12 


D 


(59) 
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Substitution  from  Equations  (59)  into  Equation  (56)  then  gives 


A  n  ( 

~  .  o  2  t ,  .  . 

*R  "  2l"j  2^^  (1  +  J 


'2  n3  / 


12  +  jo3*23 


I 


(60) 


or,  using  Equations  (34)  and  (35) , 


n„ 


W^'7 2J7-J  ”2 

3  Z1 


2tt  n^n. 


However,  for  numerical  calculations  it  is  usually  better  to  use  the 

unapproximated  version  in  Equation  (56) . 

To  obtain  'znn  and  z  n  we  use  respectively  the  expressions  for 
KU  XU 

the  0  and  X  waves  in  Equations  (57)  to  obtain  ,  B^,  Cq ,  D^;  A^,  B^, 
Cx,  D^.  We  then  have,  in  accordance  with  Equation  (56) , 


R0 


n20  B0  "JA0 
10  2tt  n10n30  D0  +  Jc0 


(62) 


Ao  n2X  BX~:!AX 


RX 


IX  2tt  alxt»3x  0x+icx 
Likewise,  using  Equation  (41),  we  have 


'10 


dz 


'IX 


‘BO 


Jq  n02(z) 


*BX 


dz 


Jq  nx  (z) 


(63) 


If  these  expressions  for  the  equivalent  complex  heights  of  reflection 
and  the  equivalent  complex  heights  of  the  bottom  of  the  ionosphere  are 
substituted  into  Equation  (16),  we  can  then  obtain  the  phase  velocity 
and  the  rate  of  attenuation  in  the  Earth-ionosphere  transmission  line 


32 


t* 


! 

< 


from  Equations  (17)  and  (18). 

There  is  little  phase-change  associated  with  intervals  of 
height  where  the  gradient  of  ionization  density  is  high,  so  that  the 
phase-change  in  Equations  (39)  is  always  numerically  small.  If  it 

is  taken  to  be  zero,  a  simple  physical  interpretation  may  be  given  for 
the  calculations  performed  in  this  section.  Figure  4  illustrates 
this  interpretation  for  an  individual  magneto-ionic  component  (0 
wave  or  X  wave)  under  day-time  conditions.  The  smallest  solution  of 
the  relevant  Equation  (2)  is  z  =  z^,  and  this  occurs  in  the  high 
gradient  on  the  under  side  of  the  D  region.  Calculation  of  reflection 
from  this  gradient  is  as  though  from  a  discontinuity  of  refractive 
index  n^-1,  where  n^  is  the  complex  refractive  index  at  z  =  z^. 
Calculation  of  reflection  from  the  high  gradient  on  the  under  side  of 
the  E  region  is  as  though  from  a  discontinuity  of  refractive  index 
n^-n^,  where  n7  and  n^  are  the  complex  refractive  indices  at  heights 
Zr,  and  z^,  the  second  and  third  solutions  of  the  relevant  Equation  (2). 
The  complex  phase-change  between  these  two  discontinuity  levels  is 

calculated  in  accordance  with  Equation  (34).  The  process  of  com¬ 
bining  the  reflections  from  both  discontinuity  levels  is  then  a 
standard  calculation  in  optics  or  in  transmission-line  theory.  However, 
it  is  necessary  to  switch  from  nearly  vertical  transmission  to  nearly 
horizontal  transmission  at  the  bottom  of  the  ionosphere;  this  occurs  in 
the  refracting  stratum  in  accordance  with  the  theory  of  GG.  It  is 
also  necessary  to  combine  the  calculations  using  the  0  and  X  wa-  es 
before  the  bounds. y  conditions  at  the  surface  of  the  Earth  are  satis- 
f  ied . 
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Chapter  3 

Figure  4.  Illustrating  simultaneous  reflection  of  the  0  or  X 
wave  from  the  D  and  E  regions  at  ELF  when  the 
relevant  Equation  (2)  has  a  triplet  of  solutions 
zv  z2,  z ^  (Z]L  <  z^,  <  Zj)  at  which  the  refractive 

indices  are  n. ,  n. ,  n,  respectively. 


.  Theory  for  five  reflection  heights 


Under  day-time  conditions  the  theory  presented  in  the  previous 
section  is  required  whenever  there  is  significant  reflection  both 
from  the  high  gradient  on  the  under  side  of  the  D  region  and  from  the 
high  gradient  on  the  under  side  of  the  E  region.  However,  the  theory 
of  the  previous  section  is  also  required  under  night-time  conditions 
unless  the  gradient  on  the  top  side  of  the  E  region  is  sufficiently 
low. 

At  night,  reflection  can  occur  both  from  the  gradient  on  the 
under  side  of  the  E  region  and  from  the  gradient  on  the  top  side  of 
the  E  region.  Moreover,  on  the  under  side  of  the  E  region  there  is 
sometimes  a  ledge  of  ionization,  somewhat  like  a  D  region  but  at  a 
greater  height.  Reflection  can  then  occur  from  the  gradient  on  the 
under  side  of  the  ledge,  from  the  gradient  on  the  under  side  of  the 
E  region,  and  from  the  gradient  on  the  top  side  of  the  E  region. 

This  behavior  manifests  itself  by  the  occurrence  of  five 
solutions  to  one  or  both  of  Equations  (2) .  Such  a  situation  is 
illustrated  in  Figure  5,  which  is  drawn  for  night-time  conditions  in  a 
middle  lattitude  using  a  model  of  the  ionosphere  of  the  type  employed 
in  the  next  chapter.  It  is  clear  that  we  need  to  be  able  to  deal 
with  five  reflection  heights. 

With  five  reflection  heights  z^,  z^,  z^,  z^  and  z,. 

(z^  <  z^  <  Zj  <  z^  <  z^),  we  employ  the  following  complex  phase-changes 
(cf  Equations  34  and  35) : 
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Illustrating  multiple  solutions  of  Equations  (2)  under 
night-time  conditions. 


fZ2 

*12-(2itAo)/  ndz, 
*23  -  (2ir/Ao)  (23  -  *2) 


»  * 


«34  -  (27rA_)  /  ndz 


(64) 


^45  *  (2^A0)(z5-z4) 


The  matrix  (38)  not  only  relates  the  values  of  E  and  H  at  height  z, 

x  y  e,  1 

to  those  at  height  z^,  but  also  those  at  height  z^  to  those  at  height 

z^  if  subscripts  1,  2  and  3  in  Equations  (39)  are  replaced  by  subscripts 

3,  4  and  5.  Moreover,  if  the  two  resulting  matrices  are  multiplied 

together,  we  have  the  matrix  that  relates  the  values  of  E  and  H  at 

x  y 

the  height  z^  to  those  at  the  height  z^.  Performance  of  this 
multiplication  is  straightforward  but  lengthy.  The  product  has  the 
form  of  the  matrix  (38)  but  with  more  complicated  expressions  for  a, 
b,  c  and  d  involving  all  the  phase-changes  appearing  in  Equations  (64). 

The  rest  of  the  calculation  in  the  previous  section  subsequent 
to  Equations  (39)  is,  however,  unchanged  up  to  Equation  (56),  which 
becomes 


~  *0  n2n4  B  -  1A 

ZR  "  21  2tt  nin3n5  D+  jC 


(65) 


where  A,  B,  C  and  D  now  have  the  values  given  in  the  Appendix  to  this 

chapter.  If  these  expressions  are  evaluated  for  the  0  and  X  waves  so 

as  to  obtain  A,.,  Bn,  C«,  D„  and  A  ,  B  ,  C  ,  D  ,  then  the  equivalent 
0  0  0  0  x  x  x  x 

complex  heights  of  reflection  for  the  0  and  X  waves  are 


n20n40 

n10n30a50 


Bo-^Ao 
Dx  + 


n2Xn4X 

nlXn3Xn5X 


V^x 

Dx+jCx 


(66) 


We  calculate  and  2^  from  these  equations,  z'gQ  and  z'BX  from 
Equations  (63),  the  eigenvalue  S  from  Equation  (16),  the  phase  velocity 
in  the  Earth-ionosphere  transmission  line  from  Equation  (17),  and  the 
rate  of  attenuation  from  Equation  (18) . 

The  phase-changes  $23  an<*  $45  associated  with  the  intervals  of 
height  where  the  gradient  of  ionization  density  is  high  are  always 
numerically  small.  If  they  are  taken  to  be  zero,  then  the  theory  in 
this  section  reduces  to  one  involving  discontinuities  of  complex 
refractive  index  at  three  levels.  The  values  of  the  discontinuities  are 
n^  -  1,  n^  -  n2  and  n,.  -  n^.  The  height  of  the  first  discontinuity  is 
z^.  The  complex  phase-change  between  the  first  discontinuity  and  the 
second  is  The  comPlex  phase-change  between  the  second  disconti¬ 

nuity  and  the  third  is  $^4*  The  effect  of  reflection  from  all  three 
discontinuities  is  calculated  as  in  optics  or  in  transmission-line 
theory . 

For  numerical  work,  however,  there  is  no  need  to  take  $23  and 
as  zero.  Moreover,  it  is  convenient  always  to  program  for  five  so¬ 
lutions  of  each  of  Equations  (2).  If,  for  one  or  both  of  the  magneto¬ 
ionic  components,  there  are  only  three  solutions  z^,  Z2  and  z^>  then 
we  put  z^  and  z^  in  the  program  equal  to  z^.  If  there  is  only  one 

solution  z_,  then  we  put  all  five  heights  in  the.  program  equal  to  z  . 

R  R 
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8.  Conclusion 


It  is  concluded  that  a  reasonable  simplified  procedure  for 
evaluating  the  effect  of  the  ionosphere  on  propagation  in  the  Earth- 
ionosphere  transmission  line  in  the  ELF  band  is  as  follows: 

(i)  Evaluate  the  heights  of  reflection  for  the  0  and  X  waves  from 
Equations  (2) . 

(ii)  At  these  heights  evaluate  the  complex  refractive  indices  n^ 
and  n^  for  vertical  propagation  of  the  0  and  X  waves. 

(iii)  At  the  same  heights  evaluate  also  the  complex  elliptical 

polarizations  Qq  and  as  defined  in  Equation  (72)  of  Chapter 
2.  Verify  that  Qq  and  Qx  are  approximately  the  same  at  all 
heights  of  reflection. 

(iv)  Calculate  the  equivalent  complex  heights  of  reflection  "z^  and 
^RX  from  E quations  (66)  using  the  expressions  for  A,  B,  C  and  D 
given  in  the  Appendix  to  this  chapter. 

(v)  Calculate  the  equivalent  complex  heights  zL,.  and  "z  of  the 
bottom  of  the  ionosphere  from  Equations  (63). 

(vi)  Substitute  into  Equation  (16)  ,  and  deduce  the  phase  velocity 
and  the  attenuation  rate  from  Equations  (17)  and  (18). 

In  the  next  chapter  this  method  is  applied  to  a  simplified 
worldwide  model  of  the  ionosphere  to  estimate  ELF  propagation  in  the 
Earth-ionosphere  transmission  line  for  a  range  of  frequencies  and  a 
range  of  ionospheric  conditions. 
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Appendix  to  Chapter  3 


When  one  of  Equations  (2)  yields  five  heights  of  reflection 

z^,  z 3,  z^  and  z^  (z^  <  z^  <  z^  <  z^  <  z^) ,  let  n^,  nj,  n^,  n^  and 

n^  be  the  complex  refractive  indices  at  these  heights  for  vertical 

propagation  of  phase  for  the  wave  concerned.  For  this  wave,  let  complex 

phase-changes  $^>  $23’  ^34  anc*  ^45  by  Equations  (64). 

The  matrix  (38)  relates  the  values  of  E  and  H'  at  height  z,  to 

x  y  1 

those  at  height  z^,  where  a,  b,  c  and  d  are  given  by  Equations  (39). 

The  same  matrix  relates  the  values  of  E  and  H  at  height  z_  to  those 

x  y  j 

at  height  z^  if  subscripts  1,  2  and  3  in  Equations  (39)  are  replaced 
by  subscripts  3,  4  and  5.  The  two  resulting  matrices  must  be  multiplied 
together. 

If  this  is  done,  the  values  of  A,  B,  C  and  D  in  Equation  (65) 
evaluate  to: 


A  ■  cos  $^  cos  $23  cos  $^4  cos  $45 


-  —  cos  $. .  cos  $„  sin  $,.  sin  $.  , 
12  23  34  45 


-  —  sin  $. ,  sin  $„,  cos  $, .  cos  $.  c 
23  34  45 


+  ~~  „  sin  $.  ,  sin  $,,,  sin  $,,  sin  $,  e 
n.n.  12  23  34  45 

2  4-  ... 

-  Oj  cos  $^2  sin  ^23  sln  ^34  cos  ^45 


-  —  cos  $, , sin  $_,  cos  .  sin  $,  e 
n^  12  23  34  45 


--sin  $12  cos  $23  sin  $34  cos  $4J 


— —  sin  $  -  cos  cos  sin  $,  c 
njO^  12  23  34  45 


(67) 
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•  iij  cos  <J>^2  c0S  VCOS  *34sin*45 
n. 

+  —  cos  $12  cos  $23  sin  <t>34  cos  <J>45 
4 

-Xin$12sin$23  cos43Asin*45 

-  ~~  sin  *12  sin  *23  sin  4>34  cos  $45 

2  4 

-  n3n5  cos  «12  sin  $23  sin  $34  sin  $43 
n,n, 

+  _i-£  cos  $i2  sLn  $23  cos  $34  cos  <i>43 
n,n, 

-  sin  4>12  cos  4>23  sin  4>34  sin  4>43 

+  *12  cos  *23  C°S  *34  C0S  *45 

2  4 

-  sin  «12  cos  $23  cos  $34  cos  4>4$ 

n- 

-  -~±-  sin  *- ,  cos  «„  sin  4>,  sin  4>43 
n3n3  a* 

+  cos  4>12  Sin  $23  cos  <?34  cos  $43 

-  cos  $12  sin  4>23  sin  4>34  sin  4>4$ 


—  sin  4,  -  sin  4’--  cos  4,.  sin  4.  , 
12  23  34  45 


n4 

+  ^  C0S  *12  <:0S  *23  Sln  *34  C0S  *45 


+  —  cos  4,  ,  cos  <J>-.  cos  4>  sin  4,  c 
n^  12  23  34  *5 


(69) 


a2n4 

D  - - —  sin  <4,  .  cos  4.,  cos  4  sin  4 

Qj  12  23  34  45 


2 

- sin  4,  .  cos  4,,  sin  4  cos  4 

n3  12  23  34  45 


—  cos  4.  .  sin  4_,  cos  4  sin  4 
n3  12  23  34  45 


-^CM*12Sln*23’l"*34C°S*M 

+  n2n4  sin  4^  sin  423  sin  4^  sin  4^ 

-  n2  sin412sin423  cos  4^  cos  4^ 

-  n4  cos  412  cos  423  sin  4^  sin  4^ 

+  cos  4,  „  cos  4  cos  4  cos  4 

12  23  34  45 


(70) 


If  these  expressions  are  evaluated  for  the  0  and  X  waves  respectively 
so  as  to  obtain  AQ,  BQf  CQ>  DQ  and  A ^  Dx,  then  the  equivalent 

complex  heights  of  reflection  for  the  0  and  X  waves  are  given  by 
Equations  (66). 

If  in  Equations  (67)  -  (70)  we  put  z 4  =  z5>  they  reduce  to 
Equations  (57).  If  in  Equations  (67)  -  (70)  we  put  zx  *  z2  =  z3  “  Z4 

z  *  z  ,  then  Equations  (66)  reduce  to  Equations  (4) • 

5  R 
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Chapter  4.  Calculations  of  ELF  Propagation  in 
Earth-ionosphere  Transmission  Line 

1.  Introduction 

In  Chapter  3  a  simplified  theory  was  derived  for  evaluating  the 
effect  of  the  ionosphere  on  ELF  propagation  in  the  Earth-ionosphere 
transmission  line.  It  is  the  objective  of  this  chapter  to  apply  the 
theory  to  a  simplified  worldwide  model  of  the  ionosphere. 

In  the  theory,  phase  (and  group)  propagation  below  the  bottom 
of  the  ionosphere  is  practically  horizontal,  and  phase  propagation  above 
the  bottom  of  the  ionosphere  (but  not,  in  general,  group  propagation) 
is  practically  vertical.  A  precise  definition  of  the  bottom  of  the 
ionosphere  has  been  given  by  Greifinger  and  Greifinger  (1978,  1979). 

It  is  a  complex  height,  the  positive  imaginary  part  of  which  is  a 
measure  of  the  absorption  experienced  in  the  lowest  layer  of  the  iono¬ 
sphere.  An  indispensable  feature  of  the  theory  presented  in  Chapter  3  is 
the  use,  in  a  high  gradient  of  refractive  index,  of  a  solution  of  the 
equations  of  propagation  of  the  type  developed  by  Greifinger  and 
Greifinger  (1978,  1979).  Use  is  also  made  in  Chapter  3  of  equivalent 
complex  heights  of  reflection  for  the  0  and  X  waves.  These  are  based 
on  the  ELF  reflection  theory  of  Booker  and  Lefeuvre  (1977)  described  in 
Chapter  2. 

The  basic  procedure  is  to  compare  the  local  vertical  gradient 
of  refractive  index  with  the  local  wavelength,  thereby  classifying 
altitude  into  intervals  where  the  gradient  is  high  and  ones  where  it  is 
low.  Where  the  gradient  is  low,  the  phase-integral  treatment  is 
adequate.  An  interval  where  the  gradient  is  high  may,  to  a  first 
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approximation,  be  replaced  by  a  discontinuity.  The  amount  of  the 
discontinuity  is  the  difference  between  the  refractive  indices  at  the 
top  and  bottom  of  the  interval  of  high  gradient,  judged  in  relation  to 
local  wavelength.  It  is  then  a  matter  of  combining  reflections 
from  the  several  discontinuities.  This  requires  calculation  of  the 
complex  phase-changes  between  the  discontinuities,  that  is,  in  the 
intervals  of  height  where  the  phase-integral  treatment  is  satisfactory. 

An  improvement  on  the  discontinuity  theory  is  obtained  by 
making  allowance  for  the  small  change  of  phase  that  occurs  in  an  in¬ 
terval  of  high  gradient.  This  may  be  done  using  the  type  of  treatment 
developed  by  Greifinger  and  Greifinger  (1977,  1979).  If  this  improve¬ 
ment  is  incorporated,  calculations  are  made  as  follows: 

(i)  Evaluate  the  heights  of  reflection  z  and  z  for  the  0  and  X 

U  X 

waves  in  accordance  with  Equations  (2)  of  Chapter  3.  These 
are  algebraic  relations  that  will  be  referred  to  as  the  reflec¬ 
tion  equations. 

(ii)  At  these  heights,  evaluate  the  complex  refractive  indices  n^ 

and  n^  for  the  0  and  X  waves  for  vertical  propagation  of  phase 
in  the  ionosphere  in  the  presence  of  the  Earth's  magnetic  field, 

(iii)  At  the  same  heights,  evaluate  also  the  complex  elliptical 

polarizations  Qq  and  for  the  0  and  X  waves.  Verify  that  Qq 
and  Qx  are  approximately  the  same  at  all  heights  of  reflection, 

(iv)  Evaluate  the  equivalent  complex  heights  of  reflection  "z^  and 
z^  for  the  0  and  X  waves  as  described  in  Chapter  3. 

(v)  Evaluate  the  equivalent  complex  heights  and  z_v  of  the 

dU  BA 

bottom  of  the  ionosphere  using  the  method  of  Greifinger  and 
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Greifinger  (1978,  1979);  see  Equations  (63)  of  Chapter  3. 

(vi)  Calculate  the  eigenvalue  S  (the  sine  of  the  complex  angle  of 

incidence  upon  the  ionosphere)  from  the  formula  (Equation  16  of 
Chapter  3) 


S 


h 


-i, 


lqo  ^bo^ro)-qx 


(1) 


(vii)  Deduce  the  ratio  of  the  velocity  c  of  light  in  free  space  to  the 
phase  velocity  v  in  the  Earth-ionosphere  transmission  line  from 
the  formula  (Equation  17  of  Chapter  3) 


c/v  »  R(S)  (2) 

and  the  rate  of  attenuation  in  nepers  per  unit  distance  from  the 
formula  (Equation  18  of  Chapter  3) 

a  -  -  (2tt/Ao)I(S)  (3) 

where  X  is  the  wavelength  in  free  space,  and  where  R  and  I 
denote  real  and  imaginary  parts. 

For  the  day-time  ionosphere,  the  method  of  calculation  takes 
into  account  reflection  from  the  gradient  on  the  under  side  of  the  D 
region  and  reflection  from  the  gradient  on  the  under  side  of  the  E 
region.  The  method  would  also  take  into  account  reflection  from  a 
gradient  on  the  top  side  of  the  E  region,  but  no  significant  such 
gradient  exists  in  the  ionospheric  model  used  in  this  chapter  for  the 
day-time  ionosphere. 

For  the  night-time  ionosphere,  the  method  of  calculation  takes 
into  account  reflection  from  the  gradient  on  the  under  side  of  the  E 
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region  and  reflection  from  the  gradient  on  the  top  side  of  the  E 
region.  However,  the  model  of  the  night-time  ionosphere  used  in  this 
chapter  incorporates  a  ledge  of  ionization  below  the  E  region,  somewhat 
like  a  D  region  but  at  a  greater  height.  Reflection  can  then  be  impor¬ 
tant  from  the  gradient  on  the  under  side  of  the  ledge  as  well  as  from 
those  on  the  under  side  and  the  top  side  of  the  E  region.  All  three 
sources  of  reflection  are  taken  into  account  in  this  chapter. 

However,  it  will  transpire  that,  at  the  low-frequency  end  of  the 
ELF  band,  a  considerable  amount  of  energy  is  transmitted  through  the 
E  region  at  night  and,  in  the  form  of  both  the  whistler  mode  and  the 
Alfvdn  mode,  is  reflected  from  the  F  region.  Reflection  of  ELF  waves 
from  the  F  region,  although  important  at  night  at  frequencies  below  the 
ionic  gyro frequency,  is  not  treated  in  this  chapter. 

When  reflection  from  the  F  region  is  important,  some  variation 
with  height  in  the  elliptical  polarizations  of  the  0  and  X  waves  occurs 
above  100  km  and  may  be  important.  However,  in  all  of  the  circumstances 
actually  treated  in  this  chapter,  we  have  found  that  the  complex 
elliptical  polarizations  of  the  magneto-ionic  components  are  practi¬ 
cally  independent  of  height  at  all  relevant  heights.  Consequently , 
Equations  (32)  of  Chapter  3  are  well  satisfied. 

For  a  magneto-ionic  component  (0  wave  or  X  wave) ,  reflecting 

strata  occur  where  1/ ( 4tt)  times  the  local  wavelength  X  is  approximately 

equal  to  the  local  scale-height  h,  defined  in  terms  of  the  profile  of 
2 

n  -  1,  where  n  is  the  ionospheric  refractive  index  for  vertical 
propagation  of  phase.  The  theory  in  Chapter  3  contemplates  a  situation 
in  which  five  such  levels  may  exist  at  the  heights  z^,  z^,  z z^  and 
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z^  (z^  '->z2  <  Z3  <  z4  <  z5^  *  These  are  the  solutions  of  the  reflec¬ 
tion  equations  (Equations  2  of  Chapter  3)  for  the  magneto-ionic  component 
concerned.  In  practice,  the  level  z^  is  in  the  gradient  on  the  under 
side  of  the  D  region  (or,  at  night,  of  a  ledge  below  the  E  region  if 
such  a  ledge  exists).  The  level  z^  is  close  to  the  maximum  of  ioniza¬ 
tion  density  in  the  D  region,  and  the  level  z^  is  in  the  gradient  on 
the  under  side  of  the  E  region.  By  day,  the  levels  z^  and  z,.  do  not 
exist  because  of  the  absence  of  a  substantial  gradient  on  the  top  side 
of  the  E  region.  By  night,  however,  the  levels  z^  and  z,.  frequently 
do  exist  and  occur  in  the  gradient  above  the  level  of  maximum 
ionization  density  in  the  E  region. 

Below  the  level  z^,  X/(4irh)  is  mostly  large  and  the  field  takes 
the  form  described  by  Greifinger  and  Greifinger  (1978,  1979).  Be¬ 
tween  the  levels  z^  and  z 2,  X/(4irh)  is  mostly  small  and  the  field  con¬ 
forms  reasonably  well  to  the  theory  for  a  slowly  varying  medium.  The 
complex  phase-change  between  these  two  levels  is 


12 


(4) 


Between  the  levels  z^  and  z^,  X/(4irh)  is  again  large  and  the  field 
takes  the  form  described  by  Greifinger  and  Greifinger,  except  that  only 
vertical  propagation  of  phase  is  now  involved.  Between  these  two 
levels  there  is  a  phase-change  $23  which  is  always  numerically  small. 
Between  the  levels  z^  and  z^  we  are  again  back  in  a  slowly  varying 
medium,  and  the  complex  phase-change  between  these  levels  is 
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4> 


34 


(2tt/X  ) 
o 


(5) 


Between  the  levels  z^  and  z^  we  return  to  a  rapidly  varying  medium 
and  the  corresponding  phase-change  is  small.  Above  the  level  z^ 

we  are  in  a  slowly  varying  medium  and  we  simply  have  a  wave  for  which 
phase  travels  vertically  upwards  indefinitely. 

If  the  small  phase-changes  $23  a°d  *45  associated  with  the 
intervals  of  high  gradient  are  taken  to  be  zero,  a  simple  physical 
interpretation  of  the  theory  in  Chapter  3  is  possible.  It  is 
illustrated  in  Figure  1  for  an  individual  magneto-ionic  component 
(0  wave  or  X  wave).  Reflection  from  the  high  gradient  on  the  under 
side  of  the  E  region  is  now  reflection  from  a  discontinuity  of 
refractive  index  n^  -  n2>  where  n2  and  n^  are  the  complex  refractive 
indices  at  the  heights  z2  and  Zy  Reflection  from  the  high  gradient 
on  the  top  side  of  the  E  region  is  reflection  from  a  discontinuity  of 
refractive  index  n^  -  n^,  where  n^  and  n^  are  the  complex  refractive 
indices  at  the  heights  z^  and  z^.  Moreover,  reflection  from  the  high 
gradient  on  the  under  side  of  the  D  region  is  reflection  from  a  dis¬ 
continuity  of  refractive  index  n^  -  1,  where  n^  is  the  complex  re¬ 
fractive  index  at  the  height  z^.  Reflection  from  the  three  dis¬ 
continuities  must  be  combined,  taking  into  account  the  complex  phase- 
change  $^2  (Equation  4)  between  the  D  level  and  lower  E  level,  and 
the  complex  phase-change  $24  (Equation  5)  between  the  lower  E  level 


and  the  upper  E  level.  This  is  a  standard  calculation  in  optics  or 
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Figure  1.  Illustrating  simultaneous  reflection  of  the  0  or  X  wave 
from  (i)  a  high  gradient  on  the  under  side  of  the  D 
region,  (ii)  a  high  gradient  on  the  under  side  of  the 
E  region,  and  (iii)  a  high  gradient  on  the  top  side  of 
the  E  region. 
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in  transmission-line  theory. 

Below  the  level  z^  in  Figure  1  it  is,  of  course,  necessary  to 
switch  from  nearly  vertical  propagation  to  nearly  horizontal  propagation 
at  the  bottom  of  the  ionosphere;  this  occurs  in  a  refracting  stratum 
as  described  by  Greifinger  and  Greifinger  (1978,  1979).  It  is  also 
necessary  to  remember  that  the  calculation  for  an  individual  magneto¬ 
ionic  component  (0  wave  or  X  wave)  does  not  satisfy  the  boundary  con¬ 
ditions  at  the  surface  of  the  Earth  because  the  field  is  elliptically 
polarized.  It  is  the  combination  of  the  0  and  X  wave  calculations  that 
is  significant  at  low  heights  and  that  satisfies  the  boundary  condi¬ 
tions  at  the  Earth's  surface. 

Figure  1  illustrates  in  a  simple  way  the  basis  of  the  calcula¬ 
tions  made  in  this  chapter.  In  the  actual  computations,  however,  the 
phase-changes  $23  an<*  ‘*’45  associated  with  the  intervals  of  high  gradient 
were  not  taken  to  be  zero.  Neither  was  any  use  made  of  approximations 
based  on  the  assumption  that  $23  anc*  ^45  are  numerically  small.  We 
simply  used  without  approximation  the  formulae  given  in  the  Appendix 
of  Chapter  3. 
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.  A  simplified  model  of  the  ionosphere 


To  apply  Equations  (1),  (2)  and  (3)  to  the  study  of  ELF  propa¬ 
gation  in  the  Earth-ionosphere  transmission  line,  we  need  a  world¬ 
wide  model  of  the  ionosphere-  In  principle,  it  is  desirable  to  use 
the  International  Reference  Ionosphere  (Rawer,  Bilitza  and  Ramakrish- 
nan  1978;  Rawer,  Ramakrishnan,  and  Bilitza  1978).  However,  we  have 
used  a  simpler  model,  but  nevertheless  one  that  is  likely  to  illustrate 
many  of  the  features  that  arise  in  worldwide  ELF  propagation. 

We  have  made  no  attempt  to  model  the  ionosphere  above  a  height 
of  150  km.  At  E  region  levels  and  below,  the  profile  of  electron 
density  that  we  have  used  under  quiet  ionospheric  conditions  is  illus¬ 
trated  in  Figure  2.  All  profiles  were  represented  as  analytic 
functions  of  height  using  the  method  of  Booker  (1977)  ;  Table  1  shows 
the  parameters  used  to  reproduce  Figure  2.  The  electron  density  is 
assumed  to  depend  on  the  sunspot  number  R  in  the  manner  shown  at  the 
head  of  Table  1.  Electron  density  was  taken  to  be  proportional  to  the 
square  root  of  the  cosine  of  the  sun’s  zenith  angle  x»  an  assumption 
that  ir.  satisfactory  for  the  E  region  but  less  so  for  the  D  region. 
Figure  2  shows  that  we  have  given  the  under  side  of  the  E  region  quite 
a  high  gradient  at  a  height  of  a  little  over  80  km.  This  was  done  in 
accordance  with  the  preliminary  version  of  the  International  Reference 
Ionosphere  (Rawer,  Ramakrishnan  and  Bilitza,  1975). 

The  profile  of  electron  density  that  we  have  used  at  night  under 
quiet  ionospheric  conditions  is  based  on  Table  2  and  is  illustrated  in 
Figure  3.  There  is  no  dependence  on  time  during  the  night,  but  there 
is  a  dependence  on  latitude  A.  Higher  values  of  electron  density  were 
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Quiet  day  profile  of  electron  density  (Figure  1  ) 
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Quiet  night  profile  of  electron  density  (Figure  2) 
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4.  Illustrating  three  degrees  of  disturbance  in  the  midday 
profile  of  electron  density  at  latitude  75  . 
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Figure  5.  Illustrating  three  degrees  of  disturbance  in  the  night' 
time  profile  of  electron  density  at  latitude  75  . 


A 

used  at  high  latitudes  because  of  auroral  effects.  Figure  3  shows  tn'at 
the  adopted  night-time  prorile  of  electron  density  has  a  high 
gradient  on  the  under  side  of  the  E  region  at  a  height  of  about  90  km, 
with  a  ledge  of  ionization  below.  The  ledge  was  added  in  the  light 
of  such  night-time  rocket  profiles  as  were  available  to  us. 

No  attempt  was  made  to  include  a  sporadic  E  layer.  Nevertheless 
such  a  layer  is  important  because  it  can  reflect  the  whistler  mode 
transmitted  into  the  ionosphere  (Barr  1977,  Pappert  and  Moler 
1978). 

In  order  to  study  the  effects  of  enhanced  ionization  below  the 

E  region  on  ELF  propagation  in  the  Earth- ionosphere  transmission  line, 

calculations  were  performed  for  the  series  of  profiles  of  electron 

density  shown  in  Figure  4.  The  profile  marked  0  is  the  profile  already 

adopted  for  latitude  75°  at  midday  under  quiet  conditions.  For  the 

profiles  marked  1,  2  and  3  in  Figure  4,  the  D  region  electron  density 

has  been  increased  by  1,  2  and  3  powers  of  ten  respectively.  A 

similar  series  of  profiles  for  night  conditions  is  shown  in  Figure  5. 

Profile  3  might  be  appropriate  to  a  polar  cap  absorption  event,  and 

profiles  1  and  2  to  events  of  this  type  having  lesser  intensity.  A 

sudden  ionospheric  disturbance  associated  with  a  solar  flare  might 

correspond  roughly  to  profile  1  in  Figure  4  extrapolated  to  the  rest  of 

the  day-time  ionosphere  on  the  assumption  that  the  electron  density  is 
U 

proportional  to  cos'x- 

At  ELF,  ions  play  a  role,  especially  at  the  lower  end  of  the 
band.  Nevertheless,  the  ion  profiles  are  not  well  established.  The 
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profiles  assumed  for  negative  ions  are  shown  in  Figure  6  (day)  and 
Figure  7  (night).  The  profiles  marked  0  are  used  under  quiet  iono¬ 
spheric  conditions.  The  profiles  marked  1,  2  and  3  are  used  in  con¬ 
junction  with  the  corresponding  electron  profiles  1,  2  and  3  in 
Figures  4  and  5.  The  positive  ion  content  is  adjusted  to  maintain 
neutrality.  It  is  assumed  that  all  ions  are  molecular  and  have  an  atomic 
mass  29.  Consequently,  the  effect  of  dissociation  on  transmission  into 
the  ionosphere  above  100  km  is  not  included. 

In  addition  to  the  profiles  of  electron  density  and  ion  density, 
knowledge  is  required  of  the  profiles  of  collisional  frequency.  The 
collisional  frequencies  of  greatest  importance  at  ELF  are  and  v^, 
of  which  the  first  is  that  for  electrons  with  neutral  molecules  and  the 
second  is  that  for  ions  with  neutral  molecules.  Their  values  in  sec 
as  functions  of  height  z  in  km  are  assumed  to  be  given  by: 


log10'Je  *  5,4  "  0,077  (2-90) 
1Og10Vi  “  4,6  ~  0,077  (z-  90) 


(6) 


At  the  levels  of  greatest  importance  for  ELF  propagation  in  the  Earth- 
ionosphere  transmission  line.  Equations  (6)  are  in  reasonable  agree¬ 
ment  with  values  given  by  Banks  and  Kockarts  (1973). 

The  Earth's  magnetic  field  is  taken  to  be  that  of  a  central 
dipole,  and  no  allowance  is  made  for  the  difference  between  geographic 
and  geomagnetic  latitude.  The  variation  with  height  in  the  magnitude 
of  the  Earth's  magnetic  field  is  taken  into  account  in  calculating  the 
the  scale-heights  h^  and  h^  defined  in  Equations  (1)  of  Chapter  3. 
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6.  Illustrating  four  day-time  profiles  of  negative  ion 
density. 
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7.  Illustrating  four  night-time  profiles  of  negative  ion 
density. 


Ill 


In  calculating  the  scale-heights,  it  is  necessary  to  evaluate 
height-derivatives  not  only  for  the  magnitude  of  the  Earth's  magnetic 
field  but,  more  importantly,  for  the  electron  density,  the  ion  densities 
and  the  collisional  frequencies.  These  derivatives  need  to  be  calculated 
analytically  from  analytic  functions  of  height.  Confusion  results  if 
unsmoothed  numerically  estimated  profiles  are  used  and  the  derivatives 
are  calculated  by  taking  differences.  Confusion  also  results  if  profiles 
are  employed  that  are  discontinuous  or  that  have  discontinuous 
derivatives.  These  difficulties  are  avoided  by  specifying  the  profiles 
using  the  method  of  Booker  (1977). 

In  calculating  the  complex  refractive  indices  and  n^  at  each 
height,  the  square  root  of  a  complex  number  has  to  be  evaluated.  One 
sign  before  the  square  root  gives  the  value  for  the  0  wave  and  the 
other  sign  that  for  the  X  wave.  The  computer  cannot  be  allowed  to 
make  this  choice  according  to  any  rules  that  happen  to  be  programmed 
into  it.  If  this  is  done,  integrals  such  as  those  in  Equations  (4) 
and  (5)  are  liable  to  use  the  0  wave  value  of  n  over  part  of  the  range 
of  integration  and  the  X  wave  value  over  the  remainder.  The  refractive 
index  for  a  particular  magneto-ionic  component  (0  wave  or  X  wave)  must 
be  calculated  in  such  a  way  that  continuity  is  maintained  over  the 
entire  range  of  heights  in  use. 

Leakage  of  energy  into  the  Earth  causes  small  additions  to  the 
values  of  c/v  and  a  calculated  in  Equations  (2)  and  (3) .  These 
additions  are  easily  made;  they  are  not  included  in  the  numerical  re¬ 
sults  summarized  in  the  remainder  of  this  chapter. 
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.  Day-time  mid-latitude  frequency  dependence 


Let  us  begin  by  considering  the  frequency  dependence  of  ELF 
propagation  at  a  middle  latitude  under  day-time  conditions.  Let  the 
sunspot  number  be  low  (R  =  10)  and  the  latitude  be  45°.  Let  the  time 
of  day  be  noon,  and  let  us  assume  equinoctial  conditions.  The  model  of 
a  quiet  ionosphere  described  in  the  previous  section  then  specifies 
all  the  profiles,  and  we  can  proceed  to  solve  Equations  (2)  of  Chapter 
3  to  obtain  the  heights  of  reflection.  The  results  are  presented  in 
Figure  8.  The  actual  points  calculated  are  connected  by  rectilinear 
segments. 

At  frequencies  below  about  30  Hz  there  is  a  single  solution  of 
the  reflection  equation  for  the  0  wave  (marked  z^q  in  Figure  8)  and 
a  single  solution  of  the  reflection  equation  for  the  X  wave  (marked 
Z3X^  *’  they  occur  che  high  gradient  on  the  under  side  of  the  E 
region.  As  the  frequency  increases,  two  additional  solutions  appear 
for  each  wave.  The  solutions  marked  z ^  a._d  z^x  occur  in  the  high 
gradient  on  the  under  side  of  the  D  region,  while  the  solutions 
marked  z^q  and  z^x  are  located  near  the  level  of  maximum  ionization 
density  for  the  D  region. 

For  each  magneto-ionic  component  the  curves  for  the  three 
heights  connect  up  to  form  a  single  continuous  curve.  These  curves 
separate  (i)  situations  in  which  the  gradient  of  ionization  density  is 
sufficiently  small  for  the  phase-integral  treatment  to  be  satisfactory 
from  (ii)  situations  in  which  the  gradient  of  ionization  density  is 
sufficiently  large  for  the  Greifinger  and  Greifinger  (1978,  1979) 
type  of  treatment  to  be  satisfactory;  below  the  refracting  stratum 
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at  the  bottom  of  the  ionosphere,  this  becomes  a  situation  corresponding 
to  a  homogeneous  isotropic  medium.  The  connections  between  the  z ^  and 
z.j  curves  in  Figure  8  in  reality  occur,  not  at  300  Hz  (the  highest 
frequency  for  which  we  have  performed  calculations),  but  in  the  VLF 
band.  The  z^  curves  in  Figure  8  drop  to  the  bottom  of  the  ionosphere 
as  the  frequency  is  increased  to  VLF  values.  In  so  doing,  they  lose 
their  reflecting  ability  (=  n  -  1  in  Figure  1).  At  Sufficiently  high 
frequencies  the  z^  level  becomes  the  level  at  which  the  limiting 
polarizations  of  downcoming  0  and  X  waves  are  determined. 

From  Figure  8  we  see  that,  at  frequencies  below  about  30  Hz, 
reflection  is  controlled  primarily  by  the  high  gradient  on  the  under 
side  of  the  E  region.  At  the  frequency  increases,  however,  reflection 
from  the  D  region  starts  to  play  a  role.  At  frequencies  above  about 
1000  Hz  reflection  from  the  D  region  is  dominant.  But  there  is  a 
substantial  range  of  frequency  in  which  reflection  from  the  high 
gradients  on  the  under  sides  of  both  the  D  and  E  regions  must  be  taken 
into  account. 

Because  the  refractive  index  is  bigger  in  the  E  region  than  in 
the  D  region,  the  discontinuity  n^  -  n^  in  Figure  1  turns  out,  in  the 
ELF  band,  to  be  about  a  power  of  ten  bigger  than  the  discontinuity 
n^  -  1.  This  means  that  the  high  gradient  on  the  under  side  of  the  E 
region  is  potentially  much  more  effective  for  blocking  transmission 
of  energy  into  the  upper  ionosphere  than  is  the  high  gradient  on  the 
under  side  of  the  D  region.  The  relative  importance  of  the  two  in 
a  particular  situation  depends  on  the  value  of  the  complex  phase-change 
shown  in  Figure  1  and  defined  in  Equation  (4). 
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The  real  and  imaginary  parts  of  are  shown  in  Figure  9 
as  functions  of  frequency  for  the  circumstances  of  Figure  8.  For  the 
X  wave,  <S^2  *-s  roughly  real  because  this  wave  is  the  propagating 
whistler  wave.  For  the  0  wave,  roughly  imaginary  and  negative 

because  this  wave  behaves  quite  evanescently .  From  Figure  9  we  see 
that  the  attenuation  of  the  whistler  wave  in  going  through  the  D  region 
has  only  reached  about  1  db  even  at  a  frequency  as  high  as  300  Hz.  The 
whistler  wave  penetrates  in  strength  to  the  E  region  in  the  VLF  band, 
and  this  is  even  truer  in  the  ELF  band.  Moreover,  Figure  9  shows  that 
the  phase-change  through  the  D  region  for  the  whistler  wave  has  only 
reached  about  30°  at  300  Hz.  For  frequencies  less  than  100  Hz,  the 
phase-change  through  the  D  region  for  the  X  wave  is  small,  and  the 
attenuation  is  negligible.  Even  for  the  largely  evanescent  0  wave 
the  attenuation  through  the  D  region  has  only  reached  about  4  db  at 
300  Hz  and  it  is  less  than  1.5  db  at  frequencies  less  than  100  Hz. 

Below  about  100  Hz  the  D  region  is  an  optically  thin,  largely  transparent 
attachment  to  the  under  side  of  the  day-time  E  region.  It  is  only  when 
the  frequency  exceeds  about  100  Hz  that  the  D  region  starts  to  play  a 
significant  role  in  the  reflection  process.  It  is  only  when  the 
frequency  exceeds  about  1000  Hz  that  the  D  region  takes  over  as  the 
dominant  reflector  in  middle  latitudes. 

We  need  to  calculate  the  effect  of  interference  between  re¬ 
flection  from  the  high  gradient  on  the  under  side  of  the  D  region  and 
reflection  from  the  high  gradient  on  the  under  side  of  the  E  region. 

As  described  in  Chapter  3,  this  involves  calculating  equivalent  complex 
heights  of  reflection  z'rq  and  for  the  0  and  X  waves.  The  values 
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Figure  9.  Illustrating  the  variation  with  frequency  of  the  complex 

phase-change  through  the  D  region  at  a  middle  latitude 
under  day-time  conditions. 
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derived  from  Figures  8  and  9  are  shown  in  Figure  10.  The  values  for 
Che  equivalent  complex  height  of  the  bottom  of  the  ionosphere  are  also 
shown.  The  quantities  z^  and  were  calculated  from  Equations  (63) of 
Chapter  3;  they  are  Indistinguishable  from  each  other  in  Figure  10,  and  are 

shown  as  the  real  and  imaginary  parts  of  a  single  complex  height  zT  . 

B 

The  quantities  z^  and  were  calcualted  from  Equations  (66)  of 
Chapter  3  using  formulae  given  in  the  Appendix  of  Chapter  3.  However, 
for  frequencies  less  than  about  100  Hz,  adequate  results  can  be  obtained 
using  the  approximate  expression  in  Equation  (61)  of  Chapter  3. 

The  equivalent  complex  height  of  reflection  includes  the  effect 
of  phase-change  at  reflection  and,  in  the  imaginary  part,  the 
effect  of  energy-loss  in  association  with  reflection,  particularly 
leakage  into  the  upper  ionosphere.  At  7.5  Hz  Figure  10  shows  that  the 
real  part  of  is  almost  identical  to  the  actual  height  of  reflection 
shown  in  Figure  8.  Reflection  of  the  whistler  wave  is  taking  place 
from  the  high  gradient  on  the  under  side  of  the  E  region  with  no 
significant  change  of  phase.  But  the  imaginary  part  of  z. D  in 
Figure  10  shows  that  there  is  significant  leakage  of  energy  into  the 
upper  ionosphere  via  the  whistler  mode.  On  the  other  hand,  the 
negligible  imaginary  part  of  z^q  in  Figure  10  shows  that,  for  the  0 
wave  at  7.5  Hz,  there  is  no  significant  leakage  of  energy  into  the 
upper  ionosphere  by  day.  This  is  because  the  0  wave  is  largely 
evanescent.  The  storage  field  above  the  level  of  reflection  in  the  E 
region  does,  however,  cause  significant  change  of  phase  at  reflection, 
and  this  results  in  the  real  part  of  z^q  in  Figure  10  being  substantially 
in  excess  of  the  actual  height  of  reflection  shown  in  Figure  8.  Except 
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Figure  10.  Illustrating,  as  functions  of  frequency  at  a  middle  latitude, 
the  complex  heights  of  reflection  and  the  complex  height 
of  the  bottom  of  the  ionosphere  under  day-time  conditions. 
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at  low  latitudes,  it  is  storage  of  0  wave  energy  above  the  level  of 
reflection  that  affects  the  phase  velocity  in  Earth- ionosphere 
transmission  line,  while  it  is  leakage  of  X  wave  energy  into  the 
upper  ionosphere  that  is  a  prime  cause  of  attenuation.  However,  some 
attenuation  in  the  Earth- ionosphere  transmission  line  is  also  caused 
by  absorption  in  the  D  region  and,  using  the  method  of  Greifinger  and 
Greifinger  (1977,  1979),  this  appears  in  Figure  10  via  the  imaginary 
part  of  the  equivalent  height  of  the  bottom  of  the  ionosphere. 

We  notice  that,  although  Figure  8  shows  the  D  region  as  coming 
into  action  as  the  frequency  increases  through  about  30  Hz  for  the 
X  wave  and  about  45  Hz  for  the  0  wave,  no  corresponding  discontinuities 
appear  in  the  equivalent  complex  heights  of  reflection  shown  in  Figure 
10.  This  is  because,  as  shown  in  Figure  9,  the  complex  phase-change 
4^2  between  the  D  and  E  region  levels  of  reflection  is  small,  and 
increases  quite  gradually  with  frequency.  The  real  part  of  'zR^ 
increases  with  frequency  up  to  300  Hz  in  Figure  10  because  of  the 
effect  of  phase-change  at  reflection. 

To  calculate  the  phase  velocity  v  and  the  attenuation  rate 
a  in  the  Earth-ionosphere  transmission  line,  the  complex  heights 
derived  in  Figure  10  have  to  be  used  in  Equation  (1),  and  hence  in 
Equations  (2)  and  (3).  Because  we  are  not  at  a  low  latitude,  the 
elliptical  polarizations  of  both  waves  are  close  to  circular  at  all 
relevant  heights,  so  that  Qq  and  in  Equations  (1)  are  almost 
exactly  -j  and  +j  respectively.  The  results  for  c/v  and  a  are  shown 
in  Figure  11. 
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Figure  11.  Illustrating  the  variation  with  frequency  of  the  velocity 
of  phase  propagation  and  the  rate  of  attenuation  at  a 
middle  latitude  under  day-time  conditions. 
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4 .  Day-time  latitude  dependence  at  the  equinox 

Let  us  now  examine  the  latitude  dependence  of  ELF  propagation 
in  the  day-time,  starting  with  a  frequency  of  45  Hz.  We  again  consider 
equinoctial  conditions.  Likewise,  we  again  let  the  time  of  day  be 
noon  and  the  sunspot  number  be  low  (R  =  10).  Solution  of  Equations 
(2)  of  Chapter  3  for  the  heights  of  reflection  yields  the  curves  shown 
in  Figure  12.  Derivation  of  the  equivalent  complex  heights  of 
reflection  and  the  equivalent  complex  height  of  the  bottom  of  the 
ionosphere  leads  to  the  results  shown  in  Figure  13. 

Figure  12  verifies  that,  for  45  Hz,  the  D  region  is  involved 
in  the  reflection  process  in  the  day-time  at  middle  latitudes.  But 
the  diagram  shows  that  the  influence  of  the  D  region  disappears  both 
at  low  latitudes  and  at  high  latitudes.  The  low  latitude  effect  is 
caused  by  the  Earth's  magnetic  field.  The  high  latitude  effect  is 
caused  by  the  decrease  in  ionization  density  at  high  solar  zenith 
angles.  Even  at  a  latitude  of  45°  the  effect  of  the  D  region  on  the 
reflection  process  is  minor  at  a  frequency  of  45  Hz  as  described  in 
Che  previous  section.  It  follows  that,  at  all  latitudes,  reflection 
at  45  Hz  in  the  day-time  is  controlled  primarily  by  the  high  gradient 
on  the  under  side  of  the  E  region. 

Figure  13  shows  that,  at  45  Hz  under  day-time  equinoctial  con¬ 
ditions,  the  real  part  of  'z_  is  nearlv  independent  of  latitude  and  has 
an  E  region  value.  The  imaginary  part  of  z'DV,  however,  increases 
steadilv  with  latitude.  This  is  because  reflection  at  the  equator  is 
metallic  in  character  but,  as  the  latitude  increases,  the  effect  of 
the  Earth's  magnetic  field  brings  into  play  leakage  of  energy  into 
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Che  upper  ionosphere  via  the  wi.istler  mode.  For  Che  same  reason 
Che  characcer  of  Che  CransmiCCed  0  wave  changes  from  mecallic  Co  evan- 
escenc  as  one  moves  away  from  Che  equaCor.  The  associated  change  of 
phase  at  reflection  causes  a  steady  increase  with  latitude  in  Che 
real  part  of  z^q. 

As  the  frequency  decreases  from  45  Hz,  the  D  region  behavior 
shown  in  Figure  12  shrinks  out  of  existence.  But  as  the  frequency 
increases  from  45  Hz,  the  D  region  behavior  spreads  in  latitude,  and 
at  150  Hz  it  extends  from  the  equator  to  polar  regions.  At  300  Hz 
the  effect  of  Che  D  region  on  the  reflection  process  is  still  moderate 
at  high  latitudes,  but  it  is  quite  substantial  at  low  latitudes.  This 
is  illustrated  in  Figure  14  which  shows,  as  functions  of  latitude  at 
300  Hz,  the  real  and  imaginary  parts  of  the  complex  phase  change 
associated  with  passage  through  the  D  region. 

At  the  equator,  vertical  propagation  at  ELF  between  the  high 
gradients  in  the  D  and  E  regions  is  metallic  in  character.  In  these 
circumstances  the  real  and  imaginary  parts  of  are  almost  exactly 

equal  in  magnitude,  and  .ne  values  are  almost  exactly  the  same  for  the 
0  and  X  waves.  This  behavior  is  seen  in  Figure  14  for  a  frequency 
of  300  Hz.  At  this  frequency  Che  attenuation  associated  with  passage 
through  the  D  region  at  the  equator  amounts  to  about  6  db.  Conse¬ 
quently  the  skin  depth  for  the  D  region  at  the  equator  at  300  Hz  is 
less  than  the  separation  between  the  levels  of  reflection  associated 
with  the  high  gradients  on  the  under  sides  of  the  D  and  E  regions. 
Reflection  from  the  D  region  is  therefore  dominant  at  the  equator  at 
300  Hz. 
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Figure  12.  Illustrating,  as  functions  of  latitude  at  45  Hz,  the 
solutions  of  the  reflection  equations  (Equations  2  of 
Chapter  3)  under  day-time  conditions. 
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Figure  23.  Illustrating,  as  functions  of  latitude  at  45  Hz,  the 
complex  heights  of  reflection  and  the  complex  height 
of  the  bottom  of  the  icr.osphere  under  day-time  conditions 
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However,  this  quickly  changes  as  one  moves  to  higher  latitudes. 
The  effect  of  the  Earth's  magnetic  field  causes  a  change  in  the  X  wave 
from  metallic  behavior  to  whistler  behavior,  and  in  the  0  wave  from 
metallic  behavior  to  evanescent  behavior.  Figure  14  shows  that,  at 
latitude  45°,  the  attenuation  of  the  evanescent  wave  in  passing  through 
the  D  region  at  300  Hz  is  down  to  4  db  and  that  for  the  whistler  wave 
is  down  to  1  db.  By  latitude  75°  these  figures  are  still  lowei ,  and 
even  the  phase-change  for  the  whistler  wave  in  going  through  the  D 
region  is  then  down  to  20°.  At  300  Hz  by  day,  reflection  is  controlled 
by  the  high  gradient  on  the  under  side  of  the  D  regions  at  low  latitudes 
but  by  the  high  gradient  on  the  under  side  of  the  E  region  at  high 
latitudes.  As  the  frequency  is  increased,  the  D  region  dominance  spreads 
from  low  latitudes  towards  the  poles;  it  extends  to  all  latitudes  at 
frequencies  above  1000  Hz. 

We  are  now  in  a  position  to  calculate  the  phase  velocity  v  and 
the  attenuation  rate  a  in  the  Earth-ionosphere  transmission  line  as 
functions  of  latitude  for  a  range  of  frequencies  in  the  ELF  band.  Re¬ 
sults  are  shown  in  Figure  15a  for  c/v  and  Figure  15b  for  a.  The 
diagrams  show  the  change  in  behavior  at  a  frequency  of  about  100  Hz 
as  the  D  region  begins  to  take  over  reflection  at  low  latitudes.  At 
300  Hz  in  Figure  15a  there  is  a  drop  in  c/v  as  the  latitude  decreases 
from  15°  to  0°.  This  is  because  the  height  of  reflection  is  falling 
to  a  D  region  value  while  the  bottom  of  the  ionosphere  is  remaining 
at  about  the  same  height;  it  is  roughly  the  square-root  of  the  ratio 
of  these  two  heights  that  controls  the  value  of  c/v  (Greifinger  and 
Greifinger  1978,  1979).  At  300  Hz  in  Figure  15b  there  is  a  continuous 
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Figure  15a.  Illustrating  the  variation  of  the  velocity  of  phase 

propagation  with  latitude  and  frequency  under  day-time 
conditions. 
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rise  in  a  as  one  moves  from  polar  regions  to  the  equator.  This  is 
caused  by  the  transition  from  E  region  reflection  to  D  region  re¬ 
flection.  At  ELF  the  discontinuity  n^  -  1  in  Figure  1  is  about  a  power 
of  ten  smaller  than  the  discontinuity  n^  -  02*  The  D  region  gradient 
lets  more  energy  through  than  the  E  region  gradient ,  and  if  this 
energy  is  absorbed  before  encountering  the  E  region,  propagation  in 
the  Earth-ionosphere  transmission  line  is  adversely  affected. 
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For  three  frequencies,  the  effect  of  increasing  the  sunspot 
numbeT  from  10  to  100  is  shown  in  Figures  16a  and  16b  for  midday  equi¬ 
noctial  conditions.  At  7,5  and  45  Hz  attenuation  in  the  Earth- 
ionosphere  transmission  line  is  slightly  reduced.  At  75  Hz  the 
attenuation  is  slightly  reduced  at  high  latitudes,  and  slightly  in¬ 
creased  at  low  latitudes.  The  last  effect  is  due  to  increased  in¬ 
fluence  of  the  D  region. 

For  two  frequencies,  the  effect  of  seasonal  variation  is  shown 
in  Figures  17a  and  17b  for  noon  conditions  with  a  low  sunspot  number. 

At  latitudes  of  75°  and  90°  under  winter  conditions  the  sun  is  below 
the  horizon  and  the  night-time  profiles  described  in  Section  2  are 
used.  A  nocturnal  profile  is  also  used  for  latitude  90°  under 
equinoctial  conditions.  Provided  that  the  sun  is  well  above  the 
horizon,  the  seasonal  variation  is  not  large.  In  comparison  with  equi¬ 
noctial  values,  attenuation  rates  are  somewhat  increased  in  the  winter 
and  somewhat  reduced  in  the  summer. 

For  a  frequency  of  75  Hz,  the  diurnal  variation  is  shown  in 
Table  3  for  equinoctial  conditions  with  a  low  sunspot  number.  The 
upper  entry  gives  c/v  and  the  lower  entry  a  in  db  Mm  \  When  the  sun 
is  low  in  the  sky,  the  attenuation  rate  in  the  Earth-ionosphere  trans¬ 
mission  line  is  increased.  However,  the  main  feature  of  diurnal 
variation  is  the  difference  between  the  day  and  night  values. 
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Figure  16a.  Illustrating  the  change  in  the  velocity  of  phase 
propagation  with  sunspot  number. 
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Figure  17a.  Illustrating  the  change  in  the  velocity  of  phase 
propagation  with  season. 
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Figure  17b.  Illustrating  the  change  in  Che  rate  of  attenuation  with 
season. 
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Diurnal  variation  of  c/v  and  a  (db  >ta  at  75  Hz 
,  (equinox,  sunspot  nuaber  =*  10) 
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6 .  Night-time  propagation  ignoring  reflection  above  100  km 

Night-time  propagation  in  the  ELF  band  presents  the  difficulty 
that  reflection  from  above  the  level  of  maximum  ionization  density  in 
the  E  region  needs  to  be  considered.  Nevertheless,  let  us  begin  by 
excluding  this  phenomenon.  We  are  then  discussing,  not  the  profiles 
of  electron  density  actually  shown  in  Figure  3,  but  ones  that  are 
modified  so  as  to  keep  the  electron  density  more-or-less  uniform  at 
heights  above  about  100  km. 

Calculations  are  then  qualitatively  similar  to  those  for  the 
day-time  ionosphere,  the  ledge  of  ionization  below  the  E  region  in 
Figure  3  playing  the  role  of  the  day-time  D  region  in  Figure  2. 

Figures  15a  and  15b  for  day-time  conditions  are  replaced  by  Figures  18a 
and  18b  for  night-time  conditions  in  circumstances  when  reflection 
from  above  the  level  of  maximum  ionization  density  in  the  E  region  is 
ignored . 

Comparison  of  Figure  18a  with  Figure  15a  shows  that  the  values 
of  c/v  are  less  by  night  than  they  are  by  day.  The  height  of  re¬ 
flection  is  greater  by  night,  but  the  height  of  the  bottom  of  the  iono¬ 
sphere  is  substantially  greater.  The  ratio  of  the  two  is  changed  so 
as  to  reduce  the  values  of  c/v. 

Comparison  of  Figure  18b  with  Figure  15b  shows  that,  at 
frequencies  below  about  100  Hz,  the  attenuation  rate  is  greater  by 
night  than  it  is  by  day.  In  both  cases  reflection  is  primarily  from 
the  high  gradient  on  the  under  side  of  the  E  region.  But  the  lower 
values  of  the  ionization  density  at  night  lead  to  lower  values  of  the 
refractive  indices,  and  consequently  to  lower  values  of  the  disconti- 
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nuities  n^-1  and  n 2_n2  Figure  !•  This  means  that  the  gradient  on 
the  under  side  of  the  night-time  E  region  is  substantially  less 
effective  in  blocking  transmission  of  energy  into  the  upper  ionosphere 
than  is  that  for  the  day-time  E  region.  In  consequence,  rates  of 
attenuation  in  the  Earth-ionosphere  transmission  line  at  frequencies 
below  about  100  Hz  are  greater  by  night  than  they  are  by  day  if  re¬ 
flection  from  the  top  side  of  the  E  region  is  ignored. 

The  same  is  true  at  frequencies  above  about  1000  Hz;  reflection 
is  then  from  the  ledge  below  the  E  region.  However,  with  the  model  of 
the  ledge  that  we  have  used,  the  transition  from  E  reflection  to  ledge 
reflection  at  night  takes  place  at  a  higher  frequency  than  does  the 
transition  from  E  reflection  to  D  reflection  by  day.  At  300  Hz  this 
results  in  the  night-time  attenuation  in  Figure  18b  being  less  than 
the  day-time  attenuation  in  Figure  15b. 

Diagrams  such  as  those  in  Figure  18a  and  18b  would  constitute 
a  reasonable  representation  of  ELF  propagation  at  night  if  this  were 
entirely  controlled  by  reflection  below  the  level  of  maximum  ionization 
density  in  the  E  region.  In  fact,  however,  this  is  not  so. 


140 


7 .  Night-time  propagation  incorporating  reflection  above  100  km 

Let  us  now  use  the  night-time  profiles  of  electron  density 
shown  in  Figure  3  without  modification  above  100  km.  The  decrease  of 
electron  density  that  we  now  have  on  the  top  side  of  the  E  region 
causes  the  reflection  equations  to  have  an  extra  pair  of  solutions 
both  for  the  0  wave  and  for  the  X  wave,  except  at  low  latitudes.  This 
produces  the  type  of  behavior  illustrated  in  Figure  1,  except  that  the 
D  region  is  replaced  by  the  ledge  on  the  under  side  of  the  E  region. 

The  degree  of  importance  of  reflection  of  a  magneto-ionic  component 
from  the  gradient  on  the  top  side  of  the  E  region  is  strongly  influenced 
by  the  complex  phase-change  4>  experienced  in  traversing  the  E 
region.  This  complex  phase-change  is  illustrated  as  a  function  of 
frequency  in  Figure  19  for  a  middle  latitude  and  a  low  sunspot  number. 

At  300  Hz  we  see  that  the  attenuation  of  the  0  wave  in  passing 
through  the  E  region  exceeds  20  db.  This  is  because  of  the  largely 
evanescent  behavior  of  the  0  wave  in  the  E  region.  In  consequence, 
reflection  of  the  0  wave  at  300  Hz  is  confined  to  the  under  side  of 
the  E  region.  However,  the  situation  is  very  different  for  the  X  wave. 
Figure  19  shows  that  attenuation  of  the  whistler  wave  in  passing  through 
the  E  region  at  300  Hz  is  less  than  0.5  db.  Hence  the  whistler  mode 
penetrates  in  great  strength  to  the  gradient  on  the  top  side  of  the 
nocturnal  E  region. 

Furthermore,  Figure  19  shows  that  the  phase-change  for  the 
whistler  mode  between  the  gradient  on  the  bottom  side  of  the  E 
region  and  the  gradient  on  the  top  side  is  close  to  180°  at  300  Hz. 

The  E  region  is  then  a  slab  with  an  optical  thickness  for  the  whistler 
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Figure  19.  Illustrating  the  variation  with  frequency  of  the  complex 

phase-change  through  the  E  region  at  a  middle  latitude 

under  night-time  conditions. 
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wave  close  to  half  a  wavelength.  This  produces  a  resonance  phenomenon. 
The  E  region  would  be  perfectly  transparent  to  the  whistler  mode 
at  the  resonant  frequency  in  the  absence  of  absorption.  Near  the 
resonant  frequency  high  field  strength  occurs  in  the  E  region  and,  in 
the  presence  of  collisional  damping,  substantial  loss  ensues. 

There  is  no  harmonic  of  this  resonant  frequency  because  re¬ 
flection  from  the  top  side  of  the  E  region  becomes  ineffective  before 
this  frequency  is  reached.  Even  for  the  fundamental  frequency,  reflec¬ 
tion  from  the  top  side  of  the  E  region  becomes  ineffective  at  low 
latitudes  because  of  the  modified  influence  of  the  Earth's  magnetic 
field.  This  is  illustrated  in  Figures  20a  and  20b,  which  present  c/v 
and  a  as  functions  of  latitude  over  the  frequency  range  30-300  Hz  for 
the  night-time  profiles  shown  in  Figure  3.  The  non-standard  behavior 
at  300  Hz  caused  by  the  resonance  phenomenon  is  immediately  apparent. 
However,  the  high  attenuation  rate  for  300  Hz  in  Figure  20b  disappears 
between  latitude  30°  and  latitude  15°.  Closer  to  the  equator  there 
is  no  significant  reflection  from  the  gradient  on  the  top  side  of  the 
E  region.  In  the  latitude  range  0°-15°  results  for  300  Hz  in  Figures 
20a  and  20b  are  the  same  as  those  in  Figures  18a  and  18b. 

Comparison  of  these  pairs  of  diagrams  shows  that  reflection 
from  the  gradient  on  the  top  side  of  the  E  region  extends  to  progres¬ 
sively  lower  latitudes  as  the  frequency  is  decreased  below  300  Hz. 
Whereas  it  disappears  between  30°  and  15°  at  300  Hz,  it  disappears  be¬ 
tween  10°  and  5°  at  150  Hz.  Moreover,  at  higher  latitudes,  inter¬ 
ference  between  the  reflections  from  the  gradients  on  the  top  and 
bottom  sides  of  the  E  region  is  non-resonant  at  150  Hz.  This  reduces 
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figure  20a.  Illustrating  the  "ariation  of  the  velocitv  of  phase 
propagation  with  latitude  and  frequency  under  night¬ 
time  conditions,  incorporating  reflection  from  che  gradient 
on  the  top  side  of  the  E  region. 
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Figure.  20b.  11’usi.rating  the  variation  of  Che  rate  of  attenuation  with 
latitude  and  frequency  'incur  nignt-c ine  conditions, 
*ncorporatinr  reflection  rrot.i  c‘ie  graoeut  on  cue  Cop 
side  of  the  E  region. 
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the  rate  of  attenuation  at  150  Hz  below  the  value  that  would  apply  in 
the  absence  of  reflection  from  the  top  side,  whereas  the  resonant  inter¬ 
ference  at  300  Hz  has  the  opposite  effect.  The  upshot  at  150  Hz  is 
a  maximum  in  the  attenuation  rate  at  a  latitude  in  the  vicinity  of 
10°  as  shown  in  Figure  20b.  At  75  Hz,  this  maximum  occurs  at  a  latitude 
in  the  vicinity  of  5°. 

Except  at  low  latitudes,  vertical  cross-sections  of  Figures  20a 
and  20b  give  the  resonance  curves  associated  with  the  E  region  resonance 
for  the  whistler  mode.  However,  because  calculations  have  not  been  made 
at  frequencies  above  300  Hz,  only  the  behavior  below  the  resonant  fre¬ 
quency  is  depicted. 

Let  us  now  turn  to  a  frequency  of  7.5  Hz,  using  the  night-time 
profiles  shown  in  Figure  3.  From  Figure  19  we  see  that  attenuation 
of  the  X  wave  in  traversing  the  E  region  at  night  is  completely 
negligible  at  7.5  Hz  (0.06  db).  Even  the  phase-change  is  only  about 
20°.  For  the  whistler  wave  at  7.5  Hz  at  night,  therefore,  the  E 
region  is  an  optically  thin,  nearly  transparent  windowpane.  Reflection 
comes  from  the  F  region.  This  is  confirmed  by  the  fact  that  the  fifth 
solution  of  the  reflection  equations  for  the  X  wave  at  7.5  Hz  at  night 
occurs  well  above  the  greatest  height  included  in  Figure  3. 

Even  for  the  0  wave  at  7.5  Hz  at  night,  we  see  from  Figure  19 
that  the  attenuation  through  the  E  region  is  only  about  3  db.  Although 
the  0  wave  is  largely  evanescent  in  the  E  region,  nevertheless  it 
tunnels  through  to  a  significant  extent  at  7.5  Hz  at  night.  Moreover, 
because  7.5  Hz  is  below  the  ionic  gyrof requency ,  the  0  wave  above  the 
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E  region  becomes  a  propagating  Alfv£n  mode.  This  too  is  reflected 
from  the  F  region,  although  significant  attenuation  is  involved.  The 
loss  of  energy  caused  by  partial  transmission  of  the  0  wave  through 
the  E  region  is  quite  significant,  and  is  a  prime  cause  of  attenuation 
in  the  Earth- ionosphere  transmission  line  at  7.5  Hz  by  night. 

Nocturnal  propagation  at  7.5  Hz  in  the  Earth-ionosphere  trans¬ 
mission  line  depends  on  reflection  from  the  F  region.  Because  we  have 
not  incorporated  an  F  region  in  our  model  of  the  ionosphere,  we  are 
unable  to  include  curves  for  7.5  Hz  in  Figure  20a  and  20b.  Curves 
for  15  Hz  are  omitted  for  the  same  reason. 

There  are  major  differences  between  (i)  the  curves  in  Figures 
18a  and  18b  and  (ii)  the  curves  in  Figures  20a  and  20b.  These 
differences  arise  from  the  assumption  in  the  latter  case  that  a 
substantial  drop  in  ionization  density  occurs  between  the  E  and  F 
regions  as  shown  in  Figure  3.  It  is  clear  that  the  magnitude  and 
other  features  of  this  drop  have  an  important  influence  on  ELF 
propagation  in  the  Earth-ionosphere  transmission  line  at  night.  They 
change  the  effect  of  the  interference  phenomenon  that  takes  place 
between  reflections  from  the  gradients  on  the  top  and  bottom  sides  of 
the  E  region.  The  extent  and  character  of  the  drop  are  not  well 
established.  The  drop  that  we  have  assumed  was  suggested  by  the 
preliminary  version  of  the  International  Reference  Ionosphere  (Rawer, 
Ramakrishnan  and  Bilitza,  1975).  The  actual  drop  may  be  different  and 
may  be  a  function  of  latitude  and  other  parameters.  Results  are 
likely  to  be  affected  by  various  phenomena  including  the  presence  of 
sporadic  E  ionization  (Barr  1977,  Pappert  and  Holer  1978). 
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8.  Ionospheric  disturbances 

Let  us  examine  the  likely  effect  on  day-time  ELF  propagation  of 
a  sudden  ionospheric  disturbance  (SID)  caused  by  a  solar  flare.  Let 
us  suppose  that  the  SID  occurs  at  midday  under  equinoctial  conditions. 

We  regard  the  SID  as  creating  a  Disturbance  1  condition  as  described  in 
Section  2.  This  increases  the  ionisation  density  by  about  a  power  of 
ten  in  the  D  region  while  leaving  it  relatively  unaffected  elsewhere. 

The  effect  on  c/v  and  a  is  shown  in  Figures  21a  and  21b  for  three 
frequencies. 

At  7.5  Hz  and  45  Hz  the  SID  lowers  the  height  of  the  bottom 
of  the  ionosphere  and  consequently  increases  the  values  of  c/v  as  shown 
in  Figure  21a.  The  same  is  true  at  75  Hz  at  high  latitudes.  At 
low  latitudes,  however,  reflection  at  75  Hz  is  transferred  from  the  E 
region  to  the  D  region  by  the  SID.  Consequently  the  ratio  of  the 
height  of  reflection  to  the  height  of  the  bottom  of  the  ionosphere  is 
reduced,  and  with  it  the  value  of  c/v.  Transfer  of  reflection  from 
the  E  region  to  the  D  region  begins  at  low  latitudes  and  spreads  to 
all  latitudes  as  the  frequency  is  increased.  For  a  quiet  ionosphere  this 
transfer  takes  place  at  a  frequency  of  about  300  Hz  (Figures  15a  and 
15b) ,  but  during  an  SID  the  transfer  of  reflection  to  the  D  region 
is  already  well  under  way  at  a  frequency  of  75  Hz. 

Figure  21b  shows  that  an  SID  increases  the  attenuation  rate  in 
the  Earth-ionosphere  transmission  line  at  7.5  Hz  only  slightly.  At 
45  Hz  the  increased  loss  is  more  substantial.  In  Figure  1,  the  en¬ 
hanced  D  region  increases  the  value  of  the  discontinuity  -  1  but 
reduces  the  value  of  the  discontinuity  n^  -  n^,  and  it  is  the  latter 
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that  is  the  more  serious  effect  at  45  Hz.  At  75  Hz  the  decrease 
in  n^  -  n^  is  also  important,  but  the  increase  in  n  -  1  makes 
itself  felt  at  low  latitudes  (Figure  21b). 

Let  us  now  consider  the  likely  effect  on  ELF  propagation  of  a 
polar  cap  absorption  (PCA)  event  caused  by  particle  emission  from  the 

sun.  At  a  latitude  of  75°  we  regard  PCA  events  of  progressively  in- 

) 

creasing  strength  as  represented  by  the  Disturbance  1,  2  and  3  con¬ 
ditions  described  in  Section  2.  The  effect  of  the  various  degrees  of 
disturbance  on  ELF  propagation  in  polar  regions  is  shown  in  Table  4  for 
three  frequencies.  At  these  frequencies  disturbance  increases 
attenuation  in  the  Earth-ionosphere  transmission  line.  The  principal 
reason  is  that,  under  quiet  ionospheric  conditions,  the  D  region  plays 
little  role  in  the  reflection  process  below  100  Hz,  and  leakage  of 
energy  into  the  upper  ionosphere  is  controlled  primarily  by  the  high 
gradient  on  the  under  side  of  the  E  region.  This  advantage  is 
weakened  under  Disturbance  1  conditions  and  destroyed  under  Dis¬ 
turbance  2  conditions.  When  the  reflection  process  has  been  firmly 
taken  over  by  ionization  below  the  E  region,  further  increase  in 
disturbance  reduces  attenuation  in  the  Earth-ionosphere  transmission 
lire.  Nevertheless,  under  Disturbance  3  conditions,  the  attenuation 
rates  at  45  and  75  Hz  are  greater  than  they  are  for  a  quiet  ionosphere — 
a  result  that  could,  however,  be  changed  by  making  the  gradient  on 
the  under  side  of  the  ionosphere  greater  for  greater  degrees  of 
disturbance.  On  the  other  hand,  at  a  frequency  of  1000  Hz,  reflection 
is  from  the  D  region  (or  at  night  from  the  ledge  on  the  under  side 
of  the  E  region)  even  under  quiet  ionospheric  conditions.  The  effect 
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RATIO  OF  VELOCITY  OF  LIGHT  TO  VELOCITY  OF  PHASE  P 


LATITUDE 

Chapter  4 

Figure  21a.  Illustrating  the  cnar.ge  in  the  velocity  of  phase 

propagation  caused  by  a  sudden  ionospheric  disturbance. 


RATE  OF  ATTENUATION  IN  DECIBEL  MEGAMETER- 


LATITUDE 

Chapter  d 

Figure  21b.  Illustrating  the  change  In  the  rate  of  attenuation  caused 
by  a  sudden  ionospheric  disturbance. 


of  disturbance  is  then  to  reduce  the  attenuation  rate  in  the  Earth- 
ionosphere  transmission  line. 


Frequency 

(hertz) 


Day  profiles  (Figs.  4  and  6) 


Night  profiles  (Figs.  5  and  7) 


(hertz) 

0 

1 

2 

3 

0 

1 

2 

3 

7.5 

' 

1.34 

0.14 

1.39 

0.17 

1.36 

0.24 

1.34 

0.20 

F  region  is 
involved 

1.34 

0.18 

1.41 

0.16 

45 

1.26 

0.43 

1.31 

0.68 

1.23 

0.84 

1.22 

0.73 

1.11 

0.39 

1.14 

0.65 

1.18 

0.86 

1.28 

0.78 

75 

1.25 

0.73 

1.29 

1.24 

1.20 

1.24 

1.20 

1.08 

1.12 

0.45 

1.14 

0.94 

1.13 

1.26 

1.25 

1.22 

Table  5 


Calculations  of  Pappert  and  Moler 
for  comparison  with  Table  4 


Frequency 

Day  profiles 

(Figs.  4  and 

6) 

Night  profiles  (Figs.  5  and  7) 

(hertz) 

0 

1 

2 

3  1 

| 

0  12  3 

7.5 

1.32 

1.35 

1.32 

1.32 

F  region  is  1.37  1.43 

0.12 

0.16 

0.21 

0.17 

involved  0.18  0.16 

45 

1.26 

1.26 

1.21 

1.21 

1.14  1.17  1.22  1.29 

0.46 

0.82 

0.73 

0.67 

0.28  0.70  0.88  0.81 

75 

1.25 

1.23 

1.18 

1.19 

1.15  1.14  1.16  1.27 

0.77 

1.47 

1.07 

1.00 

0.40  1.98  1.56  1.26 

9.  Comparison  with  the  method  of  Papper  and  Moler  (1974) 


The  method  of  calculation  described  in  Chapter  3  and  applied 
in  this  chapter  was  designed  primarily  to  reveal  the  physics  of  ELF 
propagation.  As  a  method  for  calculating  the  numerical  values  of  c/v 
and  a  it  is  only  approximate.  Let  us  compare  these  values  with  those 
obtained  by  Pappert  and  Moler  (1974)  using  a  method  in  which  the  only 
errors  should  be  those  arising  from  digitization  of  the  profiles. 

Using  their  method,  Papper  and  Moler  (private  communication)  have 
kindly  performed  calculations  corresponding  to  those  presented  in 
Table  4.  Their  results  are  shown  in  Table  3. 

Comparison  of  Tables  4  and  5  shows  that ,  while  the  numbers  are 
not  identical,  trends  are  the  same.  The  comments  about  the  effect 
of  disturbance  on  the  polar  ionosphere  made  in  the  previous  section  on 
the  basis  of  Table  4  are  unchanged  if  use  is  made  of  Table  5  instead. 

However,  there  is  a  place  in  Tables  4  and  5  where  an  appreciable 
discrepancy  occurs  in  the  rate  of  attenuation.  This  happens  for  73  Hz 
at  night  under  Disturbance  1  conditions.  In  these  circumstances,  the 
height  of  reflection  is  between  83  and  84  km,  and  this  is  on  the  sharp 
bend  of  Profile  1  in  Figure  5.  A  small  increase  in  the  estimate  of 
height  would  then  appreciably  increase  the  estimate  of  local  scale- 
height.  This  would  increase  the  rate  of  attenuation,  thereby  creating 
better  agreement  between  the  values  shown  in  Tables  4  and  5  for 
75  Hz  at  night  under  Disturbance  1  conditions.  It  should  be  remembered 
that  the  method  of  calculation  described  in  Chapter  3  assumes  that 
the  profile  of  electron  density  is  approximately  exponential  in  behavior 
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across  Che  relevent  reflecting  stratum.  This  condition  is  often 
reasonably  well  satisfied,  but  when  it  is  not,  some  error  occurs. 

A  respect  in  which  we  disagree  with  the  calculations  of  Pappert 
and  Moler  (1974)  occurs  at  low  latitudes.  They  find  a  substantial 
dependence  of  ELF  propagation  on  A,  the  angle  made  with  the  north  by 
the  horizontal  direction  of  transmission  in  the  Earth-ionosphere 
transmission  line.  They  obtain  different  results  when  A  is 
changed  by  180°,  implying  non-reciprocal  transmission.  Other  authors 
(Field  1969,  Galejs  1972)  have  also  obtained  non-reciprocal 
effects  in  the  ELF  band.  At  VLF  the  phenomenon  is  well  known 
(Budden  1961,  Galejs  1972).  Pappert  and  Moler  (1974)  encounter 
non- reciprocal  behavior  at  ELF  not  only  by  night  but  also  by  day. 

In  contrast,  we  have  obtained  no  numerically  important  variation  with 
A,  reciprocal  or  non-reciprocal,  under  any  circumstances  in  the 
ELF  band  that  we  have  investigated. 

If  non-reciprocity  exists,  it  must  be  caused  by  the  Earth's 
magnetic  field.  Moreover,  it  cannot  be  caused  at  any  level  in  the 
ionosphere  where  the  propagation  of  phase  is  vertical.  It  must  be 
caused  where  the  transition  is  taking  place  from  nearly  horizontal 
propagation  of  phase  to  nearly  vertical  propagation.  In  the  ELF 
band  (but  not  in  the  VLF  band)  this  transition  is  confined  to  a 
refracting  stratum  at  the  bottom  of  the  ionosphere.  But  in  this 
stratum  there  is  never  any  substantial  difference  between  the  re¬ 
fractive  indices  for  the  0  and  X  waves.  This  is  especially  true  by 
day;  the  height  of  the  bottom  of  the  ionosphere  is  then  illustrated 
in  Figures  10  and  13  by  the  curves  for  the  real  part  of  'z  .  By  day 
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the  refracting  stratum  occurs  where  the  gyrofrequencies  are  small 
compared  with  the  collisional  frequencies,  and  consequently  where 
the  Earth's  magnetic  field  cannot  play  an  important  role.  We  are 
therefore  unable  to  identify  a  convincing  physical  basis  for  sig¬ 
nificant  non-reciprocity  in  the  ELF  band  by  day.  By  night  some 
non-reciprocity  is  reasonable.  The  calculations  of  Pappert  and 
Moler  (1974)  do  show  less  non-reciprocity  by  day  than  by  night, 
and  perhaps  this  is  to  be  regarded  as  adequeately  consistent  with 
the  physical  argument.  But  we  doubt  it.  So  far  as  we  are  able  to 
judge,  the  computer  program  of  Pappert  and  Moler  is  giving  correct 
results  near  the  magnetic  poles  but  not  near  the  magnetic  equator. 

When  this  point  has  been  cleared  up,  the  method  of  Pappert  and 
Moler  will  give  greater  accuracy  in  the  calculation  of  c/v  and  a 
than  can  be  obtained  by  the  method  described  in  Chapter  3.  However, 
this  greater  accuracy  will  be  achieved  using  far  more  arithmetical 
operations,  and  with  considerably  less  insight  into  the  physics  of 
the  propagation  phenomena.  For  some  purposes  it  is  physical  insight 
that  is  required  rather  than  precise  arithmetic. 
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10.  The  Schumann  resonances 


As  an  example  of  the  need  for  physical  insight  rather  than  pre¬ 
cise  arithmetic,  let  us  consider  the  Schumann  resonances  associated  with 
the  Earth- ionosphere  cavity.  For  years  it  was  taken  for  granted  that 
loss  of  energy  from  the  cavity  arose  primarily  from  absorption  below 
the  E  region.  By  some  workers  this  view  still  seems  to  be  held 
(Tran  and  Polk,  1979a, b) .  Yet  the  results  presented  in  this  paper 
suggest  that  the  physics  of  the  Schumann  resonance  phenomenon  may  be 
otherwise. 

An  important  cause  for  loss  of  energy  from  the  Earth-ionosphere 
cavity  at  the  Schumann  resonant  frequency  (~7.5  Hz)  is  leakage  through 
the  nocturnal  E  region  not  merely  via  the  whistler  mode  but  also  via 
the  Alfvdn  mode.  To  recognize  that  this  is  so,  it  is  only  necessary 
to  use  Equation  (5)  to  prepare  a  diagram  of  the  type  appearing  in 
Figure  19.  High  numerical  accuracy  is  not  needed. 

Attenuation  at  the  Schumann  resonant  frequency  in  the  day 
hemisphere  is  in  all  probability  quite  low.  To  understand  how  this 
can  be,  it  is  first  necessary  to  recognize  that  reflection  occurs  in 
a  stratum  where  l/(4r)  times  the  local  wavelength  is  of  the  same  order 
of  magnitude  as  the  local  scale-height  of  the  profile.  Even  with  now 
obsolescent  profiles,  quite  a  simple  calculation  places  the  level  of 
reflection  in  the  E  region  at  the  Schumann  resonant  frequency.  With 
profiles  based  on  rocket  observations  that  exhibit  a  high  gradient  on 
the  under  side  of  the  E  region,  the  reason  for  the  low  rate  of 
attenuation  at  7.5  Hz  in  Figure  15b  can  easily  be  understood.  High 
numerical  accuracy  is  not  needed  in  order  to  ascertain  that  the 
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attenuation  rate  in  the  Earth  ionosphere  transmission  line  over  the 
day  hemisphere  can  be  quite  low  at  7.5  Hz. 

The  important  arithmetical  weaknesses  of  the  curves  presented 
in  this  chapter  for  ELF  propagation  in  the  Earth-ionosphere  transmission 
line  are  not  connected  with  the  approximations  employed  in  the  method 
described  in  Chapter  3.  They  arise  from  the  weaknesses  of  the  iono¬ 
spheric  model  to  which  the  method  is  applied  in  this  chapter.  While 
there  is  a  need  to  improve  the  model  of  the  day-time  ionosphere,  a  far 
greater  need  exists  to  improve  the  model  of  the  night-time  ionosphere. 

In  studying  the  damping  of  the  Schumann  resonances,  a  shift  in  em¬ 
phasis  is  probably  required  away  from  absorption  of  energy  in  the 
day-time  D  region  towards  leakage  of  energy  through  the  night-time 
E  region.  To  do  this,  and  to  study  other  aspects  of  ELF  propagation 
at  night,  necessitates  a  better  model  of  the  night-time  ionosphere, 
especially  above  100  km,  than  has  been  employed  in  this  chapter. 
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11.  Conclusion 


It  is  concluded  that  the  method  of  investigating  ELF  propagation 
in  the  Earth-ionosphere  transmission  line  described  in  Chapter  3  and 
applied  in  this  chapter  (i)  gives  more  physical  insight  than  any  other 
method  currently  available,  (ii)  requires  less  effort  than  any  other 
profile-based  method  currently  available,  and  (iii)  provides  reasonable 
accuracy  in  spite  of  approximations. 

However,  to  handle  reflection  from  the  F  region  at  the  low  frequency 
end  of  the  ELF  band  at  night,  further  development  of  the  method  is  needed. 
In  particular,  between  the  E  and  F  regions  some  coupling  takes 
place  between  the  0  and  X  waves  at  frequencies  less  than  the  ionic  gyro- 
frequency.  The  extent  and  consequences  of  this  coupling  need  to  be 
investigated. 
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